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The interplay of quantum fluctuations and dissipation in chains of mesoscopic superconducting 
grains is analyzed, and the results are also applied to nanowires. It is shown that in one dimensional 
arrays of resistively shunted Josephson junctions, the superconducting-normal charge relaxation 
within the grains plays an important role. At zero temperature, two superconducting phases can 
exist, depending primarily on the strength of the dissipation. In the fully superconducting phase 
(FSC), each grain acts superconducting, and the coupling to the dissipative conduction is impor- 
tant. In the SC* phase, the dissipation is irrelevant at long wavelengths. The transition between 
these two phases is driven by quantum phase slip dipoles, and is primarily local, with continuously 
varying critical exponents. In contrast, the transition from the SC* phase to the normal metallic 
phase is a Kosterlitz-Thouless transition with global character (i.e., determined by the field behavior 
at large wavelengths) . Most interesting, is the transition from the FSC phase directly to the normal 
phase: this transition, which has mixed local and global characteristics, can be one of three distinct 
types. The corresponding segments of the phase boundary come together at bicritical points. These 
behaviors are inferred from both weak and strong coupling renormalization group analyses. At in- 
termediate temperatures, near either superconductor-to-normal phase transition, there are regimes 
of super-metallic behavior, in which the resistivity first decreases gradually with decreasing tem- 
perature before eventually increasing as temperature is lowered further. The results on chains of 
Josephson junctions are extended to continuous superconducting nanowires and the subtle issue of 
whether these can exhibit an FSC phase is considered. Potential relevance to superconductor-metal 
transitions in other systems is also discussed. 



I. INTRODUCTION 



A. Motivation 



Quantum mechanical systems that are coupled to dis- 
sipative "environments" arise in many areas of physics, 
including spin dynamics in nuclear magnetic resonance 
[1, 2], damping in atomic clocks and optical interferom- 
eters [3], dephasing and decoherence in mesoscopic sys- 
tems [4] and quantum computing [5, 6]. Also in more 
conventional condensed matter contexts, dissipation has 
been argued to play crucial roles. In particular: near to 
quantum magnetic phase transitions; [7-9] in quantum 
Hall systems; [10] and in various aspects of superconduc- 
tivity, including Josephson junctions, and thin supercon- 
ducting films and wires [11-24]. The last of these is the 
primary focus of the present paper. 

For theoretical studies of dissipative effects in macro- 
scopic — and some mesoscopic — quantum systems, the 
degrees of freedom that cause the dissipation are often 
modeled as a heat bath following Caldeira and Leggett. 
[25-27] The best studied example is a single resistively- 
shunted Josephson junction (RSJJ) [28-34]. Recent ex- 
periments by Pentilla et al. [35, 36] have shown good 
agreement with the theoretical predictions. Extensions 
to arrays of RS JJs have been analyzed by several groups 
focusing on the existence and location of phase bound- 
aries between superconducting and insulating regimes of 



the set of junctions. [11, 37-46] Experimental studies of 
one and two dimensional arrays of large superconducting 
grains coupled by dissipative Josephson junctions, agree 
qualitatively with results of the theoretical analyses [46- 
57]. 

The understanding developed from studying the de- 
struction of superconductivity by quantum fluctuations 
in arrays of Josephson junctions has became a useful 
paradigm for more general considerations of quantum 
phase transitions in dissipative environments. In general, 
dissipation suppresses certain types of quantum fluctua- 
tions and thus can favor states with spontaneously broken 
symmetries, such as superconductivity. Considerable in- 
terest in further theoretical analysis of RSJJ arrays thus 
stems not only from direct experimental relevance in the 
context of superconductivity, but also from expectations 
that the concepts and approaches will be useful far more 
generally, especially for understanding universal aspects 
of quantum phase transitions that can occur at zero- 
temperature in the presence of dissipation. But to do 
this, it is crucial to take into account the small size of 
the components involved, whether they are small grains 
or individual atoms in a crystal. 

Because the models on which they are based were ini- 
tially introduced to understand macroscopic quantum 
phenomena, [25-27] most theoretical analyses of RSJJ 
arrays have assumed that the effective charges asso- 
ciated with the superconducting and normal currents 
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are perfectly mixed within the superconducting grains 
while passing separately between grains via the Joseph- 
son junctions and normal "shunts", respectively,. Such 
an approximation is reasonable for macroscopic super- 
conductors but will break down in mesoscopic or micro- 
scopic systems [58]. As we showed in Ref. [58], a con- 
sequence of the breakdown of the macroscopic paradigm 
is that the superconducting and normal fluids can effec- 
tively decouple at low energies. For the simple case of 
two junctions in series through a mesoscopic grain, we 
showed that this leads to changes in both the nature and 
location of the superconductor-normal transition that oc- 
curs. In the present paper, this analysis is extended to 
show that decoupling of the two fluids has equally impor- 
tant consequences for chains of mesoscopic grains and for 
superconducting nanowires. 

The goal of this paper is to provide a detailed anal- 
ysis of superconductor-to-normal transitions in one di- 
mensional mesoscopic systems for which dissipation plays 
a role. The main emphasis is on chains of mesoscopic 
grains that are connected both by Josephson junctions 
and some form of shunt resistance, although we also con- 
sider continuous superconducting nanowires. We will not 
discuss the possible origins of the assumed Ohmic dis- 
sipation in such systems, but rather assume that it is 
present and study its consequences. Furthermore, be- 
cause we are primarily interested in dissipative effects 
that arise simply only if the diameters of the grains or 
wires are substantially larger than atomic sizes, we will 
also ignore the limits in which the discreteness of the 
electrons or Cooper pairs becomes most essential, such 
as in Giamarchi and Schulz's treatment of superfluid-to- 
normal transitions in Luttinger liquids, [59] and phenom- 
ena associated with localization by randomness. Within 
a simple but relatively general mesoscopic model, we an- 
alyze the nature of the phases that can exist and the 
locations and character of the several types of quantum 
phase transitions that occur. In particular, we study the 
universal scaling behavior of the resistivity in the vicin- 
ity of the superconductor-to-normal transition(s). To do 
so, wc develop both strong- and weak-coupling renormal- 
ization group approaches which are tailored to deal with 
both the local resistive and the long-wavelength super- 
conducting degrees of freedom. 

B. Outline 

This paper is organized as follows. 

• Sec. I provides a general introduction. Sec. I C intro- 
duces the two-fluid approach to mesoscopic supercon- 
ducting grains, as first given in Ref. 58. Sec. ID gives 
a summary of the main results of the paper, omitting 
technical details. 

• Sec. II derives the quantum two fluid model that de- 



scribes an infinite chain of mesoscopic two-fluid grains 
as shown in Fig. 1. Sec. II A discusses the various 
possible regimes and Sec. II B gives an analysis of the 
linear electrodynamics of the model which provides in- 
tuition for the location and nature of the various tran- 
sitions. 

• In Sec. Ill we discuss the strong coupling limit of 
the chain, first deriving the quantum phase-slip rep- 
resentation of the chain (Sec. Ill A), from which we 
construct a sine-Gordon action (Sec. IIIB). Using 
the sine-Gordon action for the chain, we discuss the 
possible phases of the system (Sec. IIIC). Finally, 
in Sec. HID we derive the strong coupling RG flow 
equations for the system. These are (;onstructcd by an 
anisotropic scaling procedure, suited to the dissipative 
environment. 

• In Sec. IV wc use the RG flow equations to deter- 
mine the phase diagram of the system. Although the 
system exhibits three phases, the transitions between 
them have a variety of types. Each transition is dis- 
cussed separately in Sees. IVA-IVF. Special aspects 
of the phase diagram, such as bicritical and multicrit- 
ical points are discussed in Sees. IVG-IVI. 

• Sec. V analyzes the weak- Josephson coupling limit. 
First, we cast the action in terms of pair-tunnel events 
(Sec. VA). Then we construct the RG flow equations 
(Sec. VB), from which we obtain the weak coupling 
flow diagram (Sec. VC). 

• Sec. VI considers superconducting nanowires by con- 
sidering them as the continuum limit of the JJ chains. 

• Sec. VII presents scaling forms for the resistivity of 
the chain and discusses various interesting parameter 

regimes. 

• We conclude in Sec. VIII by reviewing the implica- 
tions of our results for various experimental systems 
and raising open questions. 

Some technical details are relegated to the Appendix. 

C. Quantum Two Fluid Description 

Mesoscopic Grains and Shunted Josephson Junctions 

The primary system we will study is shown in Fig. 
1(a): a chain of identical mesoscopic superconducting 
grains which are connected by weak links that allow the 
flow of both Cooper pairs and normal electrons. Follow- 
ing Refs. 11, 39, 40 we describe this system as a chain of 
resistively shunted Josephson junctions (Fig. 1(b)). Such 
an RSJJ chain has a natural interpretation in terms of a 
"two fluid" model: Cooper pairs that can tunnel between 
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FIG. 1; (a) A system of superconducting grains connected by weak links. The grains are coupled by flows of both supercon- 
ducting Cooper-pairs and normal electrons, (b) We model the system in (a) by a chain of superconducting grains connected 
by Josephson junctions and shunt resistors, (c) Two fluid model: each mesoscopic grain is represented as a combination of a 
superconducting and a normal parts, depicted here as separate grains. Charge relaxation between the S and TV grains is via a 
conversion resistance r. 



the superconducting grains comprise the superfluid, and 
electrons that can flow in the shunting resistors the nor- 
mal fluid. The presence of both fluids suggests consid- 
ering each grain as consisting of two (physicaUy overlap- 
ping) parts — a superconducting grain (S) and a normal 
grain (N) — as shown in Fig. 1(c). 

Changes in the super or normal charge on a grain will 
induce changes in the corresponding electrochemical po- 
tentials with the coefficients depending on both the ca- 
pacitance of the grain and the compressibilities of the nor- 
mal and superfluid components. A difference between the 
normal and the superfluid electrochemical potentials will 
lead to charge relaxation between them; we model this as 
a conversion current, /tvs, through a phenomenological 
Ohmic conversion resistance, r. 

This classical model of a chain of grains can be made 
into a quantum model straightforwardly by analogy with 
previous work, e.g., ref. [58]. The length scale is set 
by the spacing, a, between the grain centers, and vari- 
ous time scales by the normal-superfluid relaxation rate 
within a grain, the plasma frequency of the Josephson 
junctions, and the high frequency cutoffs of the supercon- 
ducting degrees of freedom (typically of order the energy 
gap) and of the dissipative processes. 



We will show that the behavior of the quantum chain 
depends crucially on the normal-superfluid conversion re- 
sistance, r. In the limit r — s- the relaxation between the 
S and AT fluids is inflnitely fast, and our model reduces to 
that of an RSJJ chain composed of macroscopic super- 
conducting grains each with a single electrochemical po- 
tential . This is the limit studied previously [42, 43, 45]. 
The opposite limit, r — > oo, describes a system with 
purely capacitative couplings between the S and the TV 
fluids. This case has been discussed in Refs. 60, 61 and a 
two dimensional version realized experimentally. [52, 62] 
The intermediate case of flnite r involves new behavior, 
which, to our knowledge, has not been analyzed previ- 
ously: it should be relevant for experiments on arrays of 
mesoscopic grains. The new physics associated with the 
interplay between the S and TV fluids dominates below 
a temperature T*, roughly proportional to the electron- 
energy-level spacing within a grain. [58] 

Nanowires 

The two fluid model can be readily generalized to a 
continuous wire that is thin enough to ignore depen- 
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dcnccs on the transverse coordinates. Classically, the cor- 
responding two-fluid model is defined by the generalized 
Josephson equation for the superfluid current, Ohm's law 
for the normal current, and a constitutive equation for 
the conversion current, 

^ = -TWVs, 

Jjv = -aVVN (1) 
^NS = I^Vn - Vs) 

with Vs and Vn the electrochemical potentials for the 
superfluid and normal electrons respectively. Note that 
for a wire, the conversion current Xns is a current per 
unit length, so 7 has dimensions conductance per unit 
length. 7, the conductivity between the normal and su- 
perfluid parts, reflects the relaxation rate for a popula- 
tion imbalance between the two fluids, as investigated, 
e.g., by Clarke et al. [63]. The current equations must 
be supplemented, as for the chain of grains, by current 
conservation laws and constitutive relations between the 
excess normal and superfluid charge densities and the 
electrochemical potentials; again these will involve com- 
pressibilities and capacitances. 

The continuum two fluid equations, Eqs. (1), can be 
made into a quantum model by defining the superfluid 
velocity as the gradient of the superconducting phase, 
imposing phase-charge conjugation, and introducing ap- 
propriate degrees of freedom to mimic the dissipation 
associated with the resistive processes. But, in addi- 
tion, the superconducting phase should be allowed to 
undergo quanf,um phase- slips (QPS): these arc implicit in 
the chain of grains but not included in the linear contin- 
uum model (for a discussion of phase-slips see Sec. III). 
To introduce phase slips, a short distance cutoff must be 
imposed, as, in contrast to the chain, there is no intrinsic 
length scale. Care must be exercised, however, when im- 
posing a cut-off on this nanowirc model. For example, a 
simple lattice regularization implicitly assumes that the 
size of quantum phase slips is also the shortest possi- 
ble distance between them. In reality the core size of a 
QPS is non-zero, but, if two phase slips occur at different 
times, their centers can be arbitrarily close to one other. 
We discuss this important subtlety in detail in Sec. VI. 

D. Overview of Results 

Before formulating the quantum model of a chain of 
grains and analyzing it in detail, we give a brief overview 
of the results obtained in this paper. 

Phases 

The J J chain has two superfluid phases, which we call 
fully superconducting (FSC) and SC* , plus a normal 



metallic phase (NOR). [The two types of superconduct- 
ing phases were discussed earlier in Refs. 42, 43, 45 and 
referred to as SC-2 and SC-1. ] In the FSC phase, which 
occurs when dissipation is strong, quantum fluctuations 
are suppressed enough that conventional superconduct- 
ing tunneling into a grain is possible, and the Josephson 
junctions behaves completely classically at zero temper- 
ature. Nevertheless the superconducting correlations are 
not truly long range, but decay algebraically with dis- 
tance. In the opposite limit, the normal phase, phase 
fluctuations are large, and number fluctuations are small, 
so that the charge of a grain is well defined and super- 
conducting correlations decay exponentially. 

In between the FSC and NOR phases is the SC* phase, 
in which both the phases and the charges of individual 
grains exhibit large fluctuations. Phase differences be- 
tween grains, however, do exhibit quasi- long range or- 
der. In this remarkable phase, tunneling into a grain 
will result in singular low energy behavior, and will be 
suppressed, which is an indication of the strong fluctu- 
ation of the phase variables on each grain in this state. 
Nevertheless a uniform supercurrent can flow through the 
chain unimpeded at low energies. The SC* phase tends 
to occur when the dissipation is weak: at low energies 
in the SC* phase, the dissipation is irrelevant with the 
quantitative values of the shunting, R, and conversion, 
r, resistances playing little role, in the determination of 
global transport properties. 

The three phase structure is somewhat analogous to 
the behavior of two-dimensional solids describable by 
disclinations: at intermediate temperatures the disclina- 
tions bind into pairs forming dislocations — the hexatic 
phase with quasi-long range orientational order — and 
at lower temperatures, these dislocations bind into pairs 
— disclination quadrupoles — forming the solid phase 
which has true orientational long range order. 

transitions 

Three types of quantum phase transitions occur be- 
tween the three phases. We denote the critical values of 
parameters for these transitions, G, L and M for global, 
local and mixed respectively. 

Separating the SC* phase and the normal phase, is 
a transition of Kosterlitz-Thouless (KT) type driven by 
unbinding of pairs of quantum phase slips. It is thus con- 
trolled by an effective space-time phase "stiffness" pa- 
rameter that we denote K, proportional to the square 
root of the ratio of the Josephson and capacitative ener- 
gies: the former favors superconductivity, while the latter 
favors localization of charge, and hence normal behavior. 
We refer to the SC* -NOR transition, as it is intrinsi- 
cally controlled by long- wavelength physics, we refer to 
as global (G). Interestingly, in this transition, the resis- 
tive shunts play almost no role: they are screened by the 
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FIG. 2: Schematic phase diagram of the quantum two fluid model of Fig. 1(c) as a function of the shunt resistance, R, and 

the quantum superconducting stiffness, K — 2tv^^J~^, for various values of the conversion resistance, r. The phase boundaries 

between the FSC, SC* , and NOR phases are of different nature: (G) indicates global (solid line), (L) local (dashed-dotted 
line), and (M) mixed (dashed lines). These come together at a multicritical point (black dot). The FSC-NOR phase boundary 
has sections with three different characters, separated by bicritical points (gray dots) Most of the phase boundaries are derived 
from the strong Josephson limit (Sec. IV), but their positions will depend on the fugacity, (, of quantum phase slips, which 
should be another axis. For small K, the weak coupling expansion (Sec. V) is needed: that portion of the FSC-NOR boundary 
is inferred from it. 



fluctuations of phase slip dipoles, i.e., a pair of phase-slip 
and anti-phase-slip occurring simultaneously on neigh- 
boring junctions. The dynamical exponent, z, that re- 
lates time or inverse-energy scales to length scales, is, for 
this transition, equal to unity: zq — I. Deviations from 
criticality are only marginally relevant, in the RG sense, 
and give rise to characteristic energy scales going to zero 
at the critical point with the exponential form typical of 
KT transitions. 

Between the two superconducting phases, the transi- 
tion is driven by dissipation. It is essentially local (L), 
being related to the superconducting-normal transition of 
a single Josephson junction. The important excitations 
that control this local transition are phase slips between 
one grain and the rest of the chain. In terms of individual 
quantum phase slips of the chain, these are equivalent to 
quantum phase slip dipoles, which, as mentioned above, 
consist of a pair of opposite sign phase slips, one on each 
side of the grain. Associated with the locality of this 
physics, there is no simple diverging length scale and the 



dynamical exponent is thus z^ — oo. How the character- 
istic energy scale goes to zero at this transition depends 
on values of resistances: the corresponding critical expo- 
nents vary continuously. 

Perhaps the most interesting transition is that which 
can occur from the FSC phase directly to the normal 
phase: this is driven by changes in the dissipation, yet 
because it also involves destruction of superconductivity, 
it has mixed (M) character with both local and global as- 
pects and involves the interplay between individual quan- 
tum phase slips, and phase slip dipoles. Related to this 
more complicated mechanism, there is more than one 
type of critical behavior — probably three types — for 
the mixed transition. 

Phase diagram 

A schematic zero-tcmperature phase diagram for the 
chain is shown in Fig. 2. It is convenient to show this as 
a function of the shunt resistance, Rq/R, and the stiff- 



6 



ncss superconducting stiffness, at various values of 
the conversion resistance, r. [Although a fourth parame- 
ter, related to QPS fugacities, is really needed as well to 
exhibit the range of possible behaviors.] Here Rq is the 
quantum resistance for Cooper pairs: 

i?Q = ^ « 6.5fcO. (2) 

For large r, r > Rq, the phase fluctuations on each grain 
arc large, and only the normal and SC* phases can ex- 
ist. In this regime, the phase diagram is simple (Fig. 
2) with the global transition between these phases at 
K — Kg ~ 4. This transition is driven by prolifera- 
tion of quantum phase slips, is analogous to the classi- 
cal Kosterlitz-Thouless transition in two dimensions, and 
has characteristic energy scale going to zero exponentially 
rapidly as the transition is approached, and a correlation 
length that diverges with the inverse of this energy scale. 

For intermediate r, Rq > r > Tc ~ jRq^ three 
phases exist, but the mixed transition between the FSC 
and normal phases is always driven by dipoles. The crit- 
ical R, Rmi varies with r and weakly with K: it is close 
to when a particular combination of R and r is -Rq . The 
energy scale goes to zero as a power of \R — Rm\ with an 
exponent, /i that varies continuously with r. The charac- 
teristic length scale on the normal side — associated with 
proliferation of individual dipoles — diverges as a power 
of i? — Rm with a continuously variable exponent that 
is not simply related to the that of the inverse energy 
scale. The FSC-SC* local transition is similarly driven 
by dipoles, and occurs when the same combination of R 
and r, is close to Rq. The energy scale similarly goes to 
zero with a continuously variable exponent. 

For r < Tc ~ \Rq - the phase diagram is far richer (Fig. 
2), and qualitatively similar to the previously studied 
r = case of infinitely fast relaxation between the nor- 
mal and superconducting electrons. The SC* -normal 
and FSC-SC* transitions are similar to the intermediate 
r regime discussed above. But the mixed-character phase 
transition between the FSC and normal phases is more 
complicated. For large-intermediate Josephson coupling, 
Kb < < 4 with Kb < 2 depending on r and R, the 
mixed transition is similar to that in the intermediate r 
regime (above). In the limit of weak Josephson coupling 
— small K — its character is different. In this regime, 
the mixed transition is most easily thought of as being 
driven — from the normal phase — by Cooper-pair tun- 
neling. It occurs when R = Rm ~ Rq — 2r with Rm 
decreasing to this value as K ^ 0. There is power law 
scaling of energy in the vicinity of this transition, with 
continuously variable exponents. 

In some regime of parameters, the small coupling and 
large-intermediate coupling phase boundaries may join 
together: if they do so, it will be in an intermediate cou- 
pling regime that is far from all three phases and thus for 
which we have no controlled methods to study. At this 



point, it is not clear whether or not these two regimes can 
join continuously. Naively , their character, particularly 
the behavior of length scales, seems rather different. But 
it is possible that they are related and then could join 
continuously. 

In the regime in which we can understand the behav- 
ior in terms of quantum phase slips, there is a section 
of the FSC-normal phase boundary which has different 
character than those discussed above. For K < Kb but 
not too small - small-intermediate - direct proliferation 
of individual QPS (rather than their proliferation caused 
by dipole proliferation) drives the transition of the the 
FSC phase. Although the behavior is more subtle and 
there are complicated crossovers, the asymptotic critical 
behavior is similar to the Kosterlitz-Thouless transition 
between the SC* and NOR phases with exponentially 
vanishing energy scale and length scale diverging as the 
inverse of this: i.e. isotropic behavior. Surprisingly, the 
location of the transition is determined by a combination 
of the low energy properties of the superconducting and 
normal — i.e. dissipative — degrees of freedom. But how 
this changes measurable asymptotic properties we have 
not worked out. 

In the low range of r in which the mixed transition can 
have more than one character, the different sections of the 
phase diagram join together at two bicritical points. The 
weak-coupling power-law regime joins up to the small- 
intermediate coupling KT-likc regime at a point that is 
neither near the K = nor the large K regime and thus 
not amenable to study by the methods we use. But the 
other bicritical point at which the dipole driven and indi- 
vidual QPS driven segments of the phase boundary come 
together can be analyzed: we discuss it briefly in Sec. 
IV G. 

The three phases come together at a multicritical point 
as shown in the flgures. The behavior near this point in- 
volves crossover from the FSC to either the SC* or normal 
phase to the left of the nearly- vertical phase boundary. 
The asymptotic critical behavior, however, is probably 
controlled by the SC* phase. We have not investigated 
this in detail. 



Renormalization group analysis 

The primary methods we use are renormalization- 
group (RG) analyses of the effective low energy action 
for the quantum two-fluid JJ chain model; both Coulomb 
gas and sine-Gordon representations of this action are 
derived and used in the various regimes. The zero- 
temperature phase diagram in various limits, and the 
nature of the quantum phase transitions more generally, 
are derived from these. The RG method naturally leads 
to detailed understanding of both the qualitative and the 
quantitative roles of various aspects of the model, for ex- 
ample, that the dissipation is irrelevant in the SC* phase 
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and near the "global" normal-to-SC* phase transition. 
For the normal-to-FSC transition, on the other hand, 
dissipation plays a key role by suppressing local fluctua- 
tions of the phase, and the characteristic length scale is 
the relaxation length between the superconducting and 
the normal fluids; as r varies from zero and infinity, this 
length-scale changes from the inter-grain spacing, a, to 
infinity, although the low energy properties of the tran- 
sition from the FSC to the SC* phase remains "local" . 

The RG flow equations can also be used to obtain the 
temperature dependence of various quantities, notably 
the resistivity close to the superconductor-to-normal 
transitions. In the vicinity of the global normal to 
SC* transition the dissipation plays little role and the 
temperature dependence of the measured total resistiv- 
ity near this quantum KT transition has been analyzed 
by other authors [64, 65] . But in the vicinity of the mixed 
normal to FSC transition, the behavior is strikingly dif- 
ferent. 

The RG analysis also gives information about 
crossovers and regimes of validity of the asymptotic 
behaviors. It is important to emphasize that over 
much of the parameter ranges, the superconductor-to- 
normal transitions are likely to be characterized by wide 
crossover regions (see Fig. 2) in which, for example, 
the resistivity can be almost temperature independent 
over extended temperature ranges. This may corre- 
spond to the "supermetallic" behavior in the vicinity of 
superconductor-to-normal transitions, that has been ob- 
served in several experiments in one dimensional [51, 66- 
69] and two dimensional systems [47, 62] . 

The RG can also be used to study finite size proper- 
ties and eff'ects of boundary conditions. For example, 
we show that when the RSJJ chain is connected to su- 
perconducting electrodes, at asymptotically low temper- 
atures in the SC* regime, the appropriate measure of the 
dissipation that controls the location of the macroscopic 
superconductor-to-normal transition is the total shunt- 
ing resistance, in contrast to the traditional picture that 
the local inter-grain shunting resistance will control the 
macroscopic behavior (see also discussion, Sec. VIII). 
[70] 

Nanowires 

Late in the paper we go briefiy from the realm of dis- 
crete grains and Josephson junctions to that of contin- 
uous superconducting wires. For these there is a subtle 
question about whether the FSC phase exists. The ap- 
proximate model we use, leads to the conclusion that it 
does not. If this is correct, then at sufficiently low tem- 
peratures, the wires will always approach the SC* or the 
normal phase, but there can be wide regimes of crossover 
in which FSC-like behavior may be observable. But as- 
pects that have been left out of the model, in particular 



aspects of charge discreteness and interference between 
quantum phase slips, may invalidate this conclusion. Pre- 
liminary indications are that these can stabilize the FSC 
phase. Its existence is thus left as an open question. 

II. MESOSCOPIC MODEL OF CHAIN OF 
GRAINS 

The system we study consists of identical mesoscopic 
superconducting grains that are connected by weak links 
which allow the flow of both Cooper pairs and normal 
electrons (Fig. 1(a)). The flow of Cooper pairs is via 
Josephson coupling between the grains and the flow of 
normal electrons via shunting resistors. As discussed 
above, it is convenient to describe this system in terms 
of a two fluid model. The superfluid is transported by 
Cooper pairs tunneling between superconducting grains, 
and the normal fluid is transported by electrons through 
the shunting resistors. Each grain is thus considered as 
a double grain with superconducting and normal parts 
as shown in Fig. 1(c). The S (superconducting) and N 
(normal) parts of each grain experience the same elec- 
trostatic potential, (p, but they will generally have differ- 
ent chemical potentials. [An analogous decoupling of the 
chemical potentials for the S and the TV fluids near phase- 
slip centers in superconducting wires in the presence of a 
transport current was discussed in Refs. 71-73.] 

The possibility of having different chemical potentials 
for the normal and superfluid components is a conse- 
quence of the mesoscopic size of the grains. The con- 
version resistance between the two fluids within a grain 
is proportional to some power of the size of a grain, and 
therefore it will provide significant dissipation only for 
small grain sizes (see discussion in Appendix A and Ref. 
[58]). The dependence of r on the size of a grain should be 
obtained from an appropriate microscopic model, which 
is beyond the scope of this paper. The sum of the 
electrostatic potential, ipi and the chemical potentials, 
MiVi! A*Si; of the i'th grain yield the total normal and su- 
perfluid electrochemical potentials, V^i and Vsi, which 
will drive the currents. 

Changes in the total charge on a grain modify its elec- 
trostatic potential via the capacitance, C, while changes 
in the normal or superfluid charges on a grain correspond- 
ingly modify the chemical potentials with the coeflicients 
the inverse compressibilities Dn ot Ds of the N and S 
fluids on an individual grain; these compressibilities thus 
have the character of "quantum capacitances" [74]. 

When the electrochemical potentials of the N and S 
fluids on a grain differ, relaxation processes will occur to 
equilibrate the two components. The simplest form for 
such relaxation is an Ohmic conversion current, InSi- 
[58] 



8 



where r is the conversion resistance within a grain. We 
assume that r remains finite even in the Umit of zero 
temperature. 

The charge relaxation time between superconducting 
and normal fluids on a grain is set by the RC time of the 
effective circuit, 

TNS = r{Ds + Dn)-^ . (4) 

Therefore our assumption of an Ohmic conversion resis- 
tance is equivalent to assuming a form for the relaxation 
rate between the normal fluid and the supcrfluid. 

The effective low-energy model of the system shown 
in Fig. 1(c) should include the charging energy for the 
grains, the Joscphson coupling energy, and appropriate 
heat bath Hamiltonians for the shunting and conversion 
resistors. We now construct the appropriate Hamiltonian 
and thence obtain the corresponding quantum action. 

We start with the charging energies. The electrochem- 
ical potentials for the superconducting and the normal 
fluids include contributions from both electrostatic and 
electrochemical "capacitances" (see Fig. 1(c)) 

Vsi = fi + DsQsi, 

VNi ^ fi + DnQni 

where ipi is the electric potential of grain i , and Qsi, Qni 
are the superconducting and the normal parts of the 
charge on grain i\ and Ds and Dm arc the inverses of 
the compressibilities of the S and N grains respectively. 
In this paper we consider a simplified model in which 
only the self-capacitance of each grain, C, is included. 
We can then write the electric potential as 

V>i = -^{Qsi + QNi)- (6) 

We expect that including mutual capacitances between 
the grains [75, 76] will not change the qualitative conclu- 
sions, although it may modify the energy scales involved. 

By integrating the electrochemical potentials in (5) 
with respect to the charge, we obtain the charging part 
of the Hamiltonian. The resulting term is 

'HQ = h{C-'+Ds)^Q%,+ 
\{C-^ + DN)T.Q%i + C-^T.QsiQNi. 



In order to write down the action for the dissipative 
and Josephson terms, we define phase-angles conjugate 
to the charges Qs and Qjv- For each grain i the supercon- 
ducting phase and the "normal phase" are defined 
via [37, 58] 

[Qjvi,Xj] = -«e^jj [Qsi,4>i] = -'^i^kj 

(8) 

[Qjvi,0j] = O [QsuXj] = ^- 

The physical interpretation of Xi follows from the obser- 
vation that its time derivative gives the electro-chemical 
potential of the normal fluid, by analogy with the Joseph- 
son relation for the superfluid [34, 37, 77] . 

The Hamiltonian term arising from the Joscphson tun- 
neling between grains can be readily written in terms of 

Hj = -EjY,cos{(t>i-<Pj), (9) 

{ij) 

where the summation is over nearest-neighbor grains, i 
and j. Ej is the Josephson coupling energy of the Joseph- 
son junctions, given by Ej = ^Ij in terms of Ij, the 
critical current of the junctions. 

The dissipative parts of the system can be described by 
heat bath Hamiltonians with appropriately chosen spec- 
tral functions. These are written as follows, 

Wdis = Wbath(i?, 2xi - 2xj) + ^bath(r, - 2xi) 

{ij) i 

(10) 

We do not give the explicit form of Wbath here, but below 
provide the corresponding effective actions obtained after 
integrating out the heat-bath degrees of freedom. The 
crucial requirement for Ohmic heat baths is that their 
density of states is linear at low frequencies. Note that 
this is the case for particle-hole excitations of a Fermi liq- 
uid, one likely source of dissipation in dirty gapless super- 
conductors especially near transitions to normal metallic 
behavior. 

Combining Eqs.(7) - (10) we construct the imaginary 
time action and partition function for the RSJJ array of 
Fig. 1(c): 



Z = J VQnVQsVc^Vxc^Y> {-S) 
S = -iiY^So drQsii'i- iY. fo dr Qm Xi + Jo dTHiQm, Qsi, (t>i,Xi) (11) 

i i 

H{QNi, Qsi, 4>ii Xi) = T^Q + + Hdis- 

I 

In the presence of Ohmic dissipation or an external cur- rent source, the phase variables (pi and Xi should be 
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periodic at t = and t = f3 with no phase twists 
by multiples of 27r allowed: in contrast to simple non- 
dissipative Hamiltonians, there is a physical distinction 
between (pi — (f)j = and 2it (see e.g. the discussion in 
Sec. IIB of Ref. 58). 

We can integrate out the {Qjvi} and {Qsi}, as these 
appear quadratically in (11), and obtain 



Z = J D(j)Dxexp{-Schain 
Schain = Sq + 5^;^ + 5^?^ + Sj 

JdrE 



Sq 



(Dn + Ds+CDnDs) 
2 I CDn l2 



l(j,i{r)-2Ur) 



5d1s = /^EEfe)¥|2x.(.„)-2x,(^n)P 
Sj = -Ej /(f dr cos (0i - (pj) 

(ij) 

(12) 

where Matsubara frequencies cOn = 2imT and Fourier 
transforms of the periodic functions (pi{T), Xii^) are used 
with f{u>n) = Jq .f(,''~)(-'''^^dT. In the limit of zero tem- 
perature, which we will primarily study, /3 E is replaced 



A. Scales, Parameters and Regimes 

There are several important energy and length scales 
in the mesoscopic model of a chain of RSJJ's, as well as 
several key dimensionless parameters. 

The superconducting energy scale is the Josephson 
coupling energy, Ej. Competing with it are (twice) the 
charging energy of a grain. 



Ec 



c 



(13) 



and the analogous non-electrostatic portions of the ener- 
gies of adding a normal electron, 



E 



DN 



N , 



or a Cooper pair. 



Eds = ^e^Ds , 



(14) 



(15) 



to a grain. In most situations of interest the energy scales 

Djv and Ds satisfy. 



Dn'-^Ds-^ 1/C , 



(16) 



so that the electrostatic energy dominates. At low ener- 
gies the D's drop out and only C is important. The ratio 



of the Josephson to the charging energy then determines 
the dimensionless quantum phase-stiffness 



K = 2-k^EjIEc 



(17) 



In the absence of dissipation, the low temperature be- 
havior is controlled by the parameter K. For large K, the 
phase differences between neighboring grains are small 
and the Josephson coupling can be approximated by 
Ej{(j)i+i —(j)if'/2. This yields the conventional quadratic 
Hamiltonian for the superconducting degrees of freedom. 
With inter-grain spacing a, the characteristic velocity of 
the phase fluctuations is 



c = 



as/ EjEc 
n ' 



(18) 



which corresponds to the Mooij-Schon velocity. [78, 79] 
In this superconducting phase, the correlations decay as 
powers of distance and imaginary time with an exponent 
proportional to K. 

The dissipation can be parametrized by the dimen- 
sionless resistances, R/Rq and v/Rq where we used the 
quantum of (Cooper pair) resistance. 



2TTn 



(19) 



these dimensionless resistances are key control parame- 
ters. When R <^ r, the effects of the two resistances — 
ii's in series and r's in parallel — becomes comparable 
at a length scale 



Aq w a^vjR 



(20) 



which plays an important role. In this small R regime, 
the dimensionless measure of the dissipation is 



Rq 



2VrR 



(21) 



When R is comparable to or larger than r, the charac- 
teristic length and dissipation measure have more com- 
plicated dependences, in particular. 



Rr 



a 



VR^ + 4ri? 



(22) 



as discussed in the next section, this can be understood 
from electro-dynamical considerations. 

The dissipative energy scale is determined by the com- 
petition between the compressibilities (which did not di- 
rectly affect the superconducting degrees of freedom) and 
the resistances. The effective resistance, R* ~ min(r, R) 
and effective capacitance. 



CsN = 



1 + CDs 



1 



Dn+Ds + CDnDs Dn + Ds 



(23) 



10 



together give the characteristic relaxation time that de- 
termines the energy scale T* parametrizing the coupling 
between the normal and superconducting currents: 



R R R R R 



R*CsN 

The energy scale, T* , can also be written in terms of the 
other energy scales. From Eq.(23), one finds that T* is 
proportional to the escape time of normal electrons from 
the grain, i.e., to the Thoulcss time of the grain. It is thus 
very large for macroscopic grains, but its existence is an 
essential property of the mesoscopic physical content of 
our model. 

At temperatures higher than T* , the dissipation across 
separate junctions is effectively decoupled. In contrast, 
for T T* , the interactions between dissipation across 
different junctions and within different grains are impor- 
tant. These interactions are crucial for the quantum dy- 
namics. At low energies, the existence of the compress- 
ibilities thus matters crucially, although the values of the 
D's do not. 

Before proceeding with more sophisticated analyses, it 
is instructive to consider the linearized dynamics in the 
presence of dissipation. As all the terms in the model ac- 
tion except the Josephson coupling are quadratic, we can 
integrate out all but the superconducting phase, </>, and 
expand the Josephson coupling energy about zero phase 
difference. This results in a rather messy form of the 
action that is given in Appendix A. In the limit of low 
frequencies and long wavelengths the dissipative effects 
are negligible and all that matters are the conventional 
superconducting parts of the action including the sup- 
pression of imaginary-time changes of the phase by the 
inverse of the total effective charging energy. These give 
rise to the simple phase modes discussed above. Never- 
theless, the fact that the short wavelength fluctuations 
arc controlled by the dissipation, makes the terms that 
appear negligible at long wavelengths also important for 
the quantum dynamics and thence the phase diagram. 

In the absence of dissipative effects, the mean square 
fluctuations of the phase difference between neighboring 
grains would be inversely proportional to the quantum 
stiffness, K. But these fluctuations are dominated by 
wavelengths of order the inter-grain spacing, and frequen- 
cies of order the Josephson junction plasma frequency 
(proportional to c/a). Thus, in actuality, the modes that 
dominate the phase fluctuations will be affected by much 
of the details of the high frequency dynamics, including 
the dissipation and the cutoff frequencies of both the dis- 
sipative and the superconducting degrees of freedom. We 
will often crudely approximate these by a high-frequency 
cutoff flo- The short wavelength processes also control 
the action of a quantum phase slip via properties of its 
core, including its space- time size and "shape"; in partic- 
ular, the "fugacity" of phase slips will be exponentially 
small when their core action is large. 





X — 


X — 


1 




X — 


X — 


(24) 




= c = 


= c 




= c = 


= c 



Phase-Slip 

FIG. 3: A phase slip creates a potential drop across the junc- 
tion, which propagates currents in the rest of the chain. 



Because of the importance of the high energy physics, 
it is not obvious what the significant dimensionless pa- 
rameters are, beyond the obvious one discussed above, 
nor what role these might play at low energies. In prac- 
tice, whether a chain of grains is in a "strong" or a "weak" 
Josephson coupling regime will be determined by many 
properties. Thus we will use these terms loosely to de- 
scribe various regimes in which the behavior simplifies: 
small QPS "fugacity", C for strong coupling, and small 
Cooper pair tunneling rates for weak coupling. Care must 
thus be exercised in c;onsi(l(;ring phase diagrams of more 
explicit models as changing one parameter can result in, 
for example, changing both the low frequency dissipation, 
and the fugacity of quantum phase slips. 

To simplify discussions of phase boundaries and the 
behavior near them, we will generally consider tuning 
the shunt resistance R and the strength of the Josephson 
coupling, either Ej itself or the QPS fugacity as a proxy 
for this. The phase diagrams in general need to be con- 
sidered as functions of both r and parameters related to 
other high energy processes as well. 



B. Circuit analysis 

In the analysis of the two junction problem in Rcf. 58 
we showed that electrical circuit properties determined 
the (lowest order) RG flows in various limits. Before 
proceeding further, we analyze the linearized electrody- 
namics of the chain model to gain insight into its behav- 
ior. 

The basis of the circuit approach is the following sim- 
ple interpretation of the dissipation-driven transition in 
a single resistively shunted Josephson junction [34, 37] . 
When the junction is in the insulating state, the Cooper 
pairs are localized on the electrodes, and the phase differ- 
ence across the junction is uncertain due to proliferation 
of quantum phase slips (QPS). In the opposite limit when 
the junction is superconducting, the phase difference be- 
tween the grains is well deflned and, phase slips will not 
occur on long time scales. A phase slip across the junc- 
tion causes a voltage burst with the Josephson relation 
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giving 

which leads to a current flow in the shunting resistor, 
Rs- It turns out that the superconducting state of the 
junction is stable only when a single QPS would cause a 
charge flow, AQ, which is more than a Cooper pair charge 
2e. From Ohm's law we then obtain AQ/2e = Rq/Rs, 
so superconductivity will be observed when Rs < Rq- 
A related analysis from the insulating state in terms of 
Cooper pair tunneling events can be used to argue for 
the same condition: Rs > Rq stabilizes the insulating 
state. 



Cooper pair tunneling 



In the normal phase the grains are only coupled resis- 

tively, although they can exhibit superconducting corre- 
lations within each grain. Thus we can consider a puta- 
tive Cooper pair tunneling from one grain to the next. 
This is very similar to the single junction case, except 
that now the total shunt resistance is R + 2r because 
of the contributions of the conversion resistances in each 
grain (sec Fig. 3). We arc interested in the effect of this 
tunneling on the SC phase difference, </> = (pi+i — (pi across 
the junction. By the Josephson relation, the change in (j) 
is given by 



A(j)= dt 



dt 



dt 



2eV{t) _ 2eV{oj = 0) 



with I the current, AQ = 2e the charge transfered, and 
Z the impedance. Thus we see that the phase difference 
between the grains will change by more than or less than 
2tt according to whether the shunting impedance, Zs = 
R+2r, is more or less than Rq. Since the transition from 
the normal phase to the FSC phase in which each grain 
has a well-defined SC phase is via this local process, it is 
not surprising that the condition for this to occur involves 



2el{w = 0)Z{u = 0) _ 2eAQZs _ ^ Zs_ 

h n ^ Rq ^ ' 
I 

this combination. As we shall see, the behavior is in fact 
more subtle due to the effects of multiple Cooper pair 
tunnels on each other. 

It is also instructive to consider the classical action 
associated with the tunneling of a Cooper pair. Since 
the phase and the Cooper pair number are canonically 
conjugate, the action will be given by 



^ I Q f , s f duj ^, . f dto I(uj) , f duj Z(ll!)\I(lo)\'^ . . 

s,. = nl ld<p ^ j dtQumt) = j -Qi.M-.) = j ^-^v^(-) = J ^^4;^ ; (27) 



note the integral of the potential energy appearing. Since 
in the limit of zero frequency the impedance is purely real 
and I{uj) — > 2e, the zero frequency part of this integral 
gives ttHZs/Rq. The finite frequency parts will be negli- 
gible if the transfer is slow; otherwise they will decrease 
the action so that 

In imaginary time, the l/{—iLj+0) becomes l/\co\ and the 
integral over frequencies diverges logarithmically at low 
frequencies or low temperatures. It is exactly the compe- 
tition between this logarithmic action and the quantum 
"entropy" — log of the range of imaginary time h/T — 
in which the event can occur — that determines whether 
the junction is superconducting, as we shall see. In real 



time, the significance of Eq. (28) is not clear, in partic- 
ular, whether ex.p{iS/h) = ex.p{in) is significant as far 
as whether or not Cooper pair tunnels can destructively 
interfere, and if they can, thereby suppressing supercon- 
ductivity for large shunt resistance. 

Effective shunting resistance in the strong Josephson 
coupling limit 

We now turn to the superconducting phase in which 
the Josephson junctions are all superconducting. Con- 
sider a phase slip across one jimction in the chain. By 
analogy with the single junction case, one would guess 
that the relevant quantity is the low frequency limit of the 
impedance, Ziuo), of the circuit parallel to it. This par- 
allel circuit involves all the other superconducting junc- 
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FIG. 4: The resistor network parallel to a single phase slip. 
The shunting resistance can be calculated using the "infinite 
resistor chain" trick. 



tions and the resistors. The total dissipation now conies 
from both the network of the Ohmic resistors and the 
Josephson junction chain itself. 

The parallel shunting impedance of a junction in the 
chain splits into two contributions: the resistance of the 
network of resistors r and R, and the impedance of the 
"telegraph line" of the chain. The first part is self ex- 
planatory; let us denote it as R^ff The second piece 
requires an explanation. When the Josephson coupling 
is strong we can assume that Josephson junctions are su- 
perconducting except for the times in which they exhibit 
a phase slip. A superconducting Josephson junction has 
the impedance of a solenoid with the ("kinetic") induc- 
tance L = h/{2e)'^ Ej^. Therefore the line of junctions 
and capacitors resembles a line of solenoids and capaci- 
tors, a simple model of a telegraph line. When a phase 
slip occurs, it creates a short-lived voltage drop on the 
junction (see Fig. 3). This pulse has two effects; the first 
is a DC current that flows through the shunt and conver- 
sion resistors, R and r and the second effect is sending 
plasma waves through the line of junctions and capaci- 
tors. The latter also behaves like a resistance as it cor- 
responds to energy being carried away by plasmons. Let 
us denote the effective resistance describing this mech- 
anism of dissipation as R^fj- The two mechanisms of 
dissipation separate at low frequencies, since the resis- 
tive contribution comes from local currents which decay 
rapidly, whereas the plasmon contributions arises from 
waves propagating at speed c. Hence, the two mecha- 
nisms are totally out of phase with each other at low 
frequencies. 

The effective resistance a phase slip feels due to the 
plasmons is 

<ff = = ■ ^ yjE^j = ^Rq (29) 

(see, for instance, Appendix 5D of Ref. 80). Note that 
R^ff include contributions from the telegraph lines on 
both sides of the given junction. 



The resistance R^ff can be found easily using the con- 
struction in Fig. 4. In the figure it is shown that the 
shunting resistance can be broken into three resistors in 
series, where the semi-infinite ladder of resistors r and R 
is replaced by an effective resistor, Rq = -^+^^'^+ ±!:Ml_ 
We thus have an effective shunting resistance of the net- 
work of Ohmic resistors 



R 



ef f — 2i?o + Ft ■ 



VR^+4:rR=^. 

a 



(30) 



The total impedance, Zs = Z{ijj = 0), shunting a junc- 
tion in the fully superconducting phase (FSC) is thus 
given by the two contributions in parallel: 



(31) 



as we shall sec, this combination controls the action of 
QPS in the FSC phase. From analogy with the single 
RSJJ case, we would expect that the FSC phase would 
become unstable to QPS when the effective inverse shunt- 
ing resistance equals Rq . As above for Cooper pair tun- 
neling, we can see that the charge transfered associated 
with the quantum phase slip of 2tt is V{uj — 0)/Zs = 
^eRqjZs- Again by analogy with the Cooper pair tun- 
neling, we can consider the action associated with the 
phase slip, finding that this is < -kTiRqIZs- Thus we 
would guess that the condition for suppression of QPS 
tunneling is when 



(32) 



However, this analogy reflects only the local physics of 
phase slips and leads to the wrong condition: neverthe- 
less, as we shall see, the correct condition involves the 
same combination. 

In Sec. VII we discuss a region of parameter space in 
the NOR phase which we call a quasi- metallic regime. In 
this region the measured resistance of the JJ chain flrst 
slowly drops as the temperature is lowered. But then, 
when a crossover temperature is reached, the resistance 
takes a sharp upturn (see Fig. 12). This effect occurs 
in the region of the NOR phase where the (wrong) local 
condition for the stability of the FSC phase to single 
phase slips (Eq. 32) is fulfllled, but the FSC is unstable 
against dipoles. 

To fully understand the effects of QPS, we also need to 
consider phase slip dipoles in which a phase slip occurs 
on one junction simultaneously with a phase slip of the 
opposite sign on another junction — say s grains away. 
In the simplest case, s = 1, this has the effect of slipping 
the phase of one grain relative to the rest of the system 
to which it is coupled. More generally, s-dipoles slip the 
phase of s consecutive grains relative to the rest of the 
chain. Because the effects of the two opposing phase slips 
cancel at long length scales, there will be no contribution 
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" f ,1 r I, r° 



: R < R 



; Ro > Rq 



Phase-Slip Dipole 

FIG. 5: Effective circuit for a phase-slip dipole with s = 1. 
Ro is the effective resistance of a semi-infinite resistor ladder 
as given above Eq. (30) . The dissipation seen by a dipole has 
no contribution from plasma waves at low frequencies. 



to the effective shunting impedance from plasmons: it 

will be entirely dissipativc. Considering all but the two 
junctions across which the slips occur, to be supercon- 
ducting and thus short circuits at low frequencies, the 
effective shunt conductance between a row of s grains 
and the rest of the chain is found to be Fig. 5, 



Rq 

R" 



2a{l - w^), 



where 



(33) 



(34) 



Of particular importance is the resistance between one 

grain and the rest of the chain that shunts the two (paral- 
lel) Josephson junctions. This is r+^{R+\/ R? + ArR) = 
Rq/13 with /3 = 2a(l - w) (see Fig. 5). 

When R <^ r. aw 1/2^/rR, (3 « Rq/r and w ~ 
1 — sjRjr so that the effective shunt resistances of strings 
of s grains only separates into two resistors in parallel — 
one at each end — for large s > Ag /a « ^JrjR. 



SC* phase and dipole proliferation 

In the FSC state individual junctions have well defined 

phase differences so that the above calculations are rele- 
vant. But in the SC* phase the superconducting phases 
of the individual grains fluctuate strongly enough to de- 
couple the superconductivity from the resistive shunts. 
In this case the energy of phase slips is primarily dissi- 
pated by the plasmons and the relevant impedance be- 
comes just the R^ff = "^Rq/K with K the dimensionless 
superconducting quantum stiffness. 

A well known result is that a Josephson junction in its 
insulating state seems like a capacitor (see, for instance. 
Appendix 5E of Ref. 80). For dipole proliferation, the 
relevant junction is the combined junction from one grain 
to the rest of the system. To show this we first note that 
the effective capacitance of a Josephson junction when 
phase slips proliferate describes the electrical response of 
the "plasma" of phase slips. This effective capacitor is 
charged whenever current tries to cross the junction itself. 




Effective Capacitance 
Due to dipoles 



FIG. 6: Phase-slip dipoles create an effective capacitance 
which inhibits current from leaving the line of Josephson junc- 
tions. This capacitance screens the normal-to-superconductor 
conversion resistance, r. But it only renormalizes the plasmar 
wave dissipation via the change of the total effective capaci- 
tance of each grain. 



For a dipole, the relevant current is the current that tries 
to leave the horizontal Josephson-j unction line in Fig. 6. 
The current leaving the junction-line into the vertically 
drawn wire is the only current that interacts with phase 
slip dipoles, hence the effective capacitor that appears 
when dipoles proliferate can get charged only by the cur- 
rent associated with Idipoies as in the Figure (6). This 
effective capacitance thus occurs between the supercon- 
ducting and normal parts of a junction, and it is in series 
to the finite conversion resistance which is important at 
low frequencies — w <C T*/h. When dipoles prolifer- 
ate, they effectively block the low-frequency conversion, 
as the effective capacitance of the dipoles dominates. 

The induced capacitance due to dipoles blocks the N-S 
conversion in each grain, and thus screens the dissipative 
interaction between phase slips: at low frequencies, this 
is essentially equivalent to setting r — > oo. But the plas- 
mon interaction is a dynamic effect. As can be seen from 
Fig. 6, the dipole-induced capacitance renormalizes the 
self capacitance C of each grain. Thus the plasmon in- 
teraction survives the proliferation of dipoles, albeit with 
a renormalizcd plasma speed. These results are formally 
derived in Sec. IV B from the effective action of the chain. 

Since dipoles disconnect the N-S conversion, in the 
SC* phase, the effective shunting resistance that is felt 
by a phase slip in the J J chain is just that from the renor- 
malizcd plasmon impedance. Therefore, 



R 



total 



K 
~2' 



(35) 



Rtotai effective local dissipation in the SC* phase. 

One would expect, in analogy to the single junction that 
when 



Rq/R, 



total 



< 1 



(36) 



the SC* phase would be unstable to QPS proliferation. 
As we shall see, the relationship between this shunting 
impedance and the stability of the SC* phase to QPS, 
is different than that for a single junction, essentially 
due to the " entropy" of the translational freedom of the 
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QPS. Nevertheless, as in the case of the FSC-NOR QPS- 
driven transition, the local condition in Eq. (36) demar- 
cates the regime with non-monotonic temperature depen- 
dence of the resistivity. In the region of the NOR phase 
where condition (36) applies, there will be a low-energy 
crossover between a decreasing resistance as temperature 
decreases, to insulating behavior (see Sec. VII). 



III. STRONG JOSEPHSON COUPLING LIMIT 

We now turn to the analysis of the low energy proper- 
ties of the chain, going well beyond the linear electrody- 
namics discussed in the previous section. We first derive 
the quantum-phase-slip Coulomb-gas representation of 
the action (Sec. Ill A and transform it into a sine-Gordon 
representation (Sec. IIIB). This enables analysis of the 
phases of the J J chain (Sec. Ill C), and an RG approach 
for the phase diagram (Sec. HID). 



A. Quantum Phase Slips and Representations of 
the Partition Function 

In the limit of strong Josephson coupling, the phase dif- 
ferences between neighboring grains will predominantly 
be localized in the vicinity of the minima of the Josephson 
potential which occur at (jji — (f)j = 27rny , with ritj an in- 
teger. Occasionally, however, the phases will depart from 
their classical superconducting form, and tunnel between 
neighboring minima with different ntj-. these events are 
quantum phase slips (QPS). Note that as discussed in 
the section on parameters, strong coupling will in prac- 
tice be defined by the rareness of QPS which depends 
on the high energy physics as well as the dimensionless 
Josephson couphng, oc Ej/Ec- 

Physically, QPS at zero temperature are caused by 
the charging energy (in Eq. 12) not commuting with 
the Josephson potential. Quantum phase slips corre- 
spond to vortex-like phase configurations in the space-r 
(imaginary-time) plane, while in real time they involve 
the launching of one dimensional plasma-waves (plas- 
mons) through the chain. 

At low temperatures, quantum phase slips may destroy 
superconducting coherence in the JJ array, and their po- 
tential role in this way makes them the basic excita- 
tions in terms of which the low energy physics can most 
readily be described when the local superconductivity is 
"strong" . We can thus try to expand in the QPS fugacity 
in this strong Josephson coupling limit. As we shall see, 
one also has to consider phase slip dipoles: bound pairs 
of QPS of opposite signs. 



Coulomb Gas Representation 

To analyze the effects of QPS in the partition function, 
it is convenient to use a Villain transformation. As usual 
first writing 

exp {/ rfrEj [cos {(t>i+i{T) - 0i(r)) - 1]} 

« E ex.p\-JdT^[(t>i+i{T)-ct>i{T) + 2niii{T)f} 

(37) 

The Villain transformation breaks the cosine function of 
phase differences down to a sum over its troughs with 
/ii(T) an integer valued function labeling the troughs in 
the Josephson potential of the junction between the i + 
I'th and the i'th grains. A phase slip on the j'th junction 
corresponds to a sudden change of fij by ±1, thus we can 
write it in terms of a density of discrete QPS: 

PQPSU, t) = = ^Pni5j,j^6{T - Tm), (38) 

m 

with the m'th QPS with "charge" Pm = ±1 occurring at 
imaginary time Tm and on junction jm- Periodic bound- 
ary conditions on the original phases, (/)i(T = 0) = (/)i(r = 
(3) imply an integrated neutrality condition for each junc- 
tion: 

/ PQPsij,r)dT = Vj. (39) 

Using Eq. (37) in the partition function, Eq. (12), 
makes the action quadratic in (^j as well as in the other 
fields, all of which can be integrated out to obtain the 
partition function solely in terms of the QPS configura- 
tions. The action of a set of N QPS with integer charges 
{Pm} at space time locations {{xm = (ijm,Tm)} , with 
(integer) jm labeling junctions, has the form: 

1 ^ 

Sn {{Pm}) = 2 PmPnG{Xm - Xn, - T„), (40) 

where G{x, r) is the interaction between phase slips. The 
QPS partition function is then: 

Zqps ^JZ^"" exp [Sn {{pm})] , (41) 

{Pm} 

with ( the fugacity of QPS; has units of frequency: in 

the absence of interactions between them, it would be 
the rate at which phase slips of each sign occur across a 
single junction. Its bare value is obtained by considering 
the action of the instanton which describes the short time 
motion of the phase variables following a phase slip in jji . 
The sum in Eq. (41) is over distinguishable QPS config- 
urations with the restriction that for each junction the 
total QPS charge is zero - a neutrality condition. Note 



15 



that strictly speaking, because of this local neutrality, 
with m and n on different junctions G{xm — Xn,Tm — Tn) 
is only the finite part of the interaction: the infinite parts 
that would arise if the local neutrality condition were not 
satisfied have been subtracted out. 

The full form of the interaction between the QPS is 
very complicated and has several different regimes. As 
discussed in the introduction, the temperature (energy) 
scale that arises from the dissipation and the finite size 
of the grains, T*, plays an important role. At tempera- 
tures higher than T* dissipation on separate junctions is 
effectively decoupled. By contrast, at low temperatures, 
T < T* , there is considerable interaction between QPS 
on neighboring junctions; we focus here and henceforth 
on this low temperature regime. 



Low Temperature Limit 

In the low temperature limit, T <^ T*, the effective 
interactions between QPS is given by 



G{x, t) « i^log 



To 



+ ae-l"l/^«log^, (42) 



where the quantum stiffness K = 2-K^fEjjEc-, (as- 
suming Dm,scC <C 1); the Mooij-Schon velocity c = 
(j^ -\/EjEc ^ g^jjjj ^jjg strength of the dissipative interactions, 



Rq 



(43) 



were all introduced earlier, and the characteristic length 
scale is 



log(l/u;) 



with 



R - V4;-/? + i?2 
2r 



(44) 



(45) 



which is approximately 1 — -y/ R/r when i? <C r. The 
origin of Aq is the decay length of currents in the resistive 
network shown in Fig. 4, and can also be inferred from 
the continuum equations (1). 

The discreteness of the chain is important even at large 
length scales, particularly when Ag is comparable to a. 
The strength of the dissipative interaction arises from 



^ / 

|fc|<fc„ 



dk 



2n r {2 -2 cos ka) + R ' 



(46) 



This yields simply the inverse of the shunt resistance for 
a single junction, as found in Sec. II B. [In the con- 
tinuum approximation with a sharp momentum space 
cutoff at n/a the QPS interaction would yield a = 



arctan (tTi,/^) : this is close to Eq.(43) over the whole 
range but does not correctly correspond to the shunt re- 
sistance. Similarly, in the continuum approximation, Aq 
would always have its asymptotic form, a\/r/R.] 

The cutoff for the plasma interaction is given by 
To ~ a/c = h/\/EjEc- The cutoff for the dissipative 
interaction in the strong coupling limit, on the other 
hand, is tq ~ min{hRQ/{R + 2r)Ej, h/y/EjEc}- When 
R + 2r ^ Rq the cutoff time for the local interaction is 
lower than that for the isotropic plasmon-related inter- 
action. This may give rise to additional crossovers in the 
temperature range REj/Rq > T > y/EjEc- However, 
we will only analyze the Coulomb-gas action, Eq. (42), 
for energy scales and temperatures < ^jEjjEc- 

The two logarithmic interactions between QPS in Eq. 
(42) have very different physical origins. The first part, 
which is isotropic in space-time, is present for l-fl dimen- 
sional X-Y models even in the absence of dissipation; it 
gives rise to a quantum Kosterlitz-Thouless transition, as 
originally discussed by Bradley and Doniach [75] (see also 
Ref. [81]). The second term in Eq. (42) originates from 
the Ohmic dissipation. It is logarithmic in the time sepa- 
ration of phase-slips. Only the limit r = (Ag = 0) has, 
to our knowledge, been analyzed previously [42, 43, 45]. 
In this limiting case the dissipative interaction of Eq. 
(42) takes the form: 



Glocal{x, t) 



a<5(x)llog^. 



(47) 



More generally, for non-zero r, the length scale, Ag, is the 
range of normal currents that are induced when a phase 
slip occurs and thus, effectively, the dissipative size of 
a QPS: Ag, controls the exponential fall-off with spatial 
separation of the dissipative interactions between QPS's. 
Not surprisingly, Ag plays a particularly important role 
in the physics of finite length chains [70]. 



B. Sine-Gordon Representation 

From the Coulomb gas partition fimction, Eq. (40), 
we derive the sine-Gordon representation of the strong 
Josephson coupling limit. The Coulomb gas representa- 
tion provides an expression for the probabilities of specific 
configurations of phase slips, whereas the sine-Gordon 
representation is more amenable to a renormalization 
group analysis. We introduce two separate Hubbard- 
Stratonovich fields, 9j and ipj located on the junctions 
between the grains + 1), in order to decouple the 
two contributions to the interactions between QPS of Eq. 
(42) (for details see Appendix 5 A. 2 of Ref. 80): in terms 
of the physical variables, ip decouples the normal degrees 
of freedom, and 9 + tp the superconducting degrees of 
freedom. This transformation yields the dual action of 
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the model of Eqs. (12): 

5duai = Jo'' dr 2^ - Oi? + !^ dr ^ E(S.^i)' 

i i 

-ClodTj:cos{e, + ^Pi) 

(48) 
(49) 

(50) 



The effective stiffnesses are 



K = K/c-- 



and 



J = Kc 



aEc' 

aEj 
2Trh ' 



Phase-slip Dipoles 

In our recent work on pairs of shunted Josephson junc- 
tions [82], we demonstrated the important role of QPS 
dipoles: these are instantons consisting of two opposite 
sign QPS that occur almost simultaneously on nearby 
junctions. The simplest phase-slip dipole, on neighboring 
junctions, is thus an event in which the phase of a single 
grain slips by 27r relative to the chain. In non-dissipative 
XY models such dipoles disappear when short wavelength 
fluctuations arc integrated out isotropically in space and 
time, and the closely-spaced dipoles act only to renormal- 
ize the stiffness k and thereby the interaction between 
other QPS. 

For granular systems with dissipation, the physics on 
the scale of the grains plays an essential role, and a bet- 
ter approach is to perform the RG procedure in time 
only. This means that dipoles remain as independent 
degrees of freedom and should be considered explicitly 
along with individual QPS. [Although it is clumsy to do 
so, we will show later that the results for the isotropic 
Kosterlitz Thouless transition can be recovered from this 
anisotropic RG. ] Dipoles can be included in the sine- 
Gordon representation by adding the following term to 
the action 

'^dipole ~ 

oo 

EVsEJdT cos m (r) + (r) ) - {6^+, (r) + (r) )] 

5=1 i 

= E E / dT7]s cos [D,9i{T) + D,^i{T)] 
i s 

(51) 

with r]s the dipole fugacity for a dipole with separation sa 

between its positive and negative QPS; we will refer to ,s, 
loosely, as the moment of the dipole; is the difference 
operator of distance s: DgOi = 6i+s — ^i- 

The argument of the cosine function in Eq. (51) places 
two phase slips of opposite signs on junctions a distance 



sa apart. Otherwise, this term is completely analogous to 
the cosine term in Eq. (48). Note that the dipole term re- 
sembles the Josephson term in the original action for the 
J J array, Eq. (12), but it is written in the dual variables. 
The similarity is explained by the following description 
of the effect of the two terms: cos{^i — removes a 
Cooper-pair at site j, and creates a Cooper pair at site i. 
Similarly, cos [{Oi{T) + ^i{T)) — {djir) + V'j ('''))]) creates 
a phase-slip on bond i, and an anti-phase-slip on bond j. 

Initially, there arc no phase slip dipoles, and rjs = for 
all s. However, upon coarse-graining in time, dipoles will 
form and have to be analyzed as independent entities, on 
the same footing as single phase slips. 



C. Phases of the JJ Array 

The three phases of the JJ array can be described in 
terms of phase slips. The simplest is the normal phase 
(NOR), in which phase slips proliferate, and phase co- 
herence is lost. In the SC* phase, phase-slips and anti- 
phase-slips bind into (neutral) dipoles and the chain be- 
comes globally superconducting although the phase fluc- 
tuations on each grain are large enough that the local 
phase is no longer well defined and the dipoles still pro- 
liferate. The phase fluctuations are strong enough that 
there is no quasi long-range order of the phase. Never- 
theless uniform supercurrent can flow unimpeded in the 
chain. 

In the fully superconducting phase (FSC), the phase 
slip dipoles bind into quadrupolcs and annihilate at low 
energy scales: the phase difl'erences on each junctions are 
then well defined, there is quasi long-range order, and 
the chain is locally, as well as globally, superconducting. 
We first study these phases and the phase transitions 
between them via a renormalization group analysis of 
the sine-Gordon representation of the QPS action. 

The two superconducting phases were discussed be- 
fore, and were dubbed respectively [42, 43, 45, 46]. The 
phase diagram of the RSJJ chain for r = was discussed 
by Chakravarty et al. [28] using the self-consistent har- 
monic approximation to treat the model of Eq. (12). 
While this analysis is appealing due to its simplicity, it 
misses an important distinction between the two distinct 
superconducting phases. A contemporary publication by 
Korshimov [42, 43], analyzed the Dyson equation of the 
sine-Gordon representation Eq.(48) and found the two 
superconducting phases, dubbing them SC-1 and SC-2 
. Unfortunately, this approach was not very transpar- 
ent and does not give a simple physical picture of the 
superconductor-to-normal transition. Fazio et al. [46] 
obtained a similar phase diagram to Korshunov via nu- 
merical calculations. In all these papers only the r = 
case was considered. In this paper we analyze the model 
of Eq. (12) for general r and provide further insight into 
the nature of the superconductor-to-normal transitions. 
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For the case r = 0, we reproduce a phase diagram derived 
by Korshunov in Refs. 42, 43. We should emphasize that 
our RG analysis is the first to allow an investigation of 
both phase boundaries and critical phenomena in the en- 
tire phase diagram (with the exception of intermediate 
Josephson strength, where, to date, no analytical method 
applies). Having the RG framework allows, for the first 
time, to discuss both finite size and finite temperature 
effects of the dissipative J J chain. 



D. Renormalization Group Flows for Strong 
Josephson Coupling 

As we showed above, in the limit of strong Josephson 
coupling the RSJJ chain of Fig. 1(c) can be described by 
a spatially discrete sine-Gordon model 



Z = I £) [6*] [■)/'] e~'^'*"'''f^''*l~'^'"p°''''^''*l. 



(52) 



We now derive and analyze the RG flow equations in 
this strong .Josephson coupling limit, treating the QPS 
fugacities as small, but at the same time allowing any 
values of the shunt and conversion resistances. 

The flows of the anisotropic RG we use, are produced in 
the sine-Gordon representation by integrating over high 
frequency modes, 9^ and ■)/''', in a frequency shell, f2 — 
dO. <\uj\ <ri, of 0i(r) and ipiir), and then rescaling the 



time as t 



1 



T SO that O returns to its 



original value, and the differential flow parameter is 



dn 



(53) 



But, crucially, lengths are not rescaled and the spatial 
discreteness is kept. In doing this, however, there ap- 
pears a characteristic length scale, A(f2), that plays an 
important role because of the anisotropic RG procedure: 
in terms of the (renormalized) phase velocity, c, which 
grows under renormalization because of the rescaling of 
time this is 



ail 



(54) 



Appendix 5.B of Ref. 80 describes in detail the deriva- 
tion of the RG-flow equations perturbatively in powers 

of the fugacities; we will need these up to second order. 

Integrating out the fast modes, effectively increases the 
short time cut-off scale from tq to To+dro = To{l+dfl/fl) 
so that part of the phase slip action is absorbed into 
its core yielding a linear suppression of the phase-slip 
fugacities. But there is also an increase in the fugacities 
from rescaling time as they represent the rate of phase 
slips. Combining these two effects wc get: 



1 



dn ( 
n \ 



1 



K 



where K 



[1 - 2a(l - w') - K{1 - e-'"/^)] } 

(55) 



The effective dissipation for the s-dipoles is 
2a(l - w^) , 



(56) 



which can be seen to be the inverse of the shunt resistance 
for the pair of (here parallel) Josephson junctions that 
couple a row of s grains to the rest of the chain, as derived 
in Sec. II B. As it is necessary to treat the smallest 
dipoles, s = 1, specially, we define 



(57) 



The K dependent term in the renormalization of rjs re- 
quires some explanation. For dipoles with large s, the two 
constituent QPS will be separately renormalized by the 
phase fluctuations that have been integrated out: thus 
the factor proportional to K in this limit. In terms of 
the circuit analysis, this dissipation can be thought of as 
arising from the frequency plasmons that propagate 
from the two opposite phase slips before annihilating far 
away. In contrast, dipoles with small s will couple only 
weakly to the low energy phase modes and thus K has lit- 
tle effect on their renormalization. The energy-dependent 
length scale A separates these two regimes. 

In principle, there should also be energy-scale depen- 
dence of the renormalization of but this can be ignored 
as long as A is substantially larger than a. We will as- 
sume this henceforth at the cost of ignoring only initial 
transient renormalizations (which will in any case only be 
part of the other relatively high energy scale processes we 
have ignored). 

Nonlinear contributions to the renormalization of C 
and the {r]s} come from second order processes. One 
such process is the combining of two phase-slips of oppo- 
site signs into a dipole. The two phase-slips have to occur 
on at almost the same time (within tq of each-other) on 
Josephson-junctions separated by s grains. This process 
gives a renormalization of ijs 

Vs^Vs+ C'^ro {2aw^ + i^e— Z^) ; (58) 

again note the s dependence of the K term: the dipoles 
are only formed with appreciable amplitude when the 
spacing between them is less than A; for convenience we 
define the factor 



u{e) = e 



(59) 



by analogy with w. 

A dipole can combine with a phase-slip on one of the 
junctions of the dipole to form a phase slip on the other 
junction. This adds second order contributions to for 
each s: 

C ^ C + CVs^To [2a{l - w') + K{1 - u')] . (60) 

Two overlapping dipoles can combine to make a single 
dipole with moment the sum of the constituent moments. 
This yields 
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dn 1 ^ 



V\a\V\s-a\ 



2a(l + w' - wl'^l - wl^-'^l) + K{1 + u' 



,|s-<rh 



(61) 



r 



where the sum runs over all a to take into the account 

both relative orientations of the dipolcs. 

Another important contribution comes from combin- 
ing two phase-slips of opposite signs on the same pair 
of Joscphson junctions with a small delay between them 
[1/At > n(l — dQ/Q)]. Their polarization gives rise to a 
renormalization of the parameter J, which is related to 
the superconducting stiffness of the Joscphson junctions. 
This yields a contribution to the quadratic action: 



V dr 



(62) 



In addition, the rescaling of imaginary time — but not 

space — leads to a rescaling of the parameters, k 
K (I — dil/n), and J ^ J {I + d^l/^) which is reflected 
in the renormalization of the velocity c and hence A oc c; 
at linear order in the fugacities, K = \fKJ is not renor- 
malized. 



The differential RG flow equations for J, C, and {r^s} 
can now be obtained by combining the terms from Eqs. 
(55, 58, 60, 61, 62) 



dl ~ 



C f 1 - - a) + 2C E VsToa{l - w'), 



^=r?, [l-2a(l-«;^)] + CVo 



(63) 



2aw* -\- \/ kJ 



+ 5 E '?<t?7|s-<7|2Q!(1 ^-w^ -w 



s _ ,„|o-| _ y,|s-o-| 



kl 



r 



valid in the limit of small fugacities C, and r\s and for low 

energy scales. As we shall see, for some purposes these 
flow equations are not sufHcient even in this limit: the 
dependence of the renormalizations of and by dipoles on 
their moments relative to A(i?) — whose renormalization 
is determined since A oc s/TJk — need to be included 
even though, formally, these disappear in the low energy 
limit of interest. 

For other purposes, we need to consider the case of 
large \r]^ but small C,. 



IV. PHASE DIAGRAM AND PHASE 
TRANSITIONS FOR STRONG JOSEPHSON 
COUPLING 

Equation (48), the action for phase slips, allows a com- 
plete description of the phase diagram in the strong cou- 
pling limit and controlled expansions in the vicinity of 
the phase transitions for some ranges of parameters. 

When the fugacity of individual QPS, C,, is relevant, 
the superconductivity should break down completely and 
the chain exhibit normal behavior. As ^ grows, it will 
also induce phase slip dipoles. If ^ is irrelevant, on the 



other hand, the chain should — at least naively — be su- 
perconducting, but which superconducting phase obtains 
depends on the behavior of the dipoles. 

The simplest phase in the strong coupling limit is the 
fully superconducting phase (FSC). In this phase, the 
fugacities {/y,,} and C are all irrelevant about the Gaus- 
sian fixed line of the sine-Gordon model, Eq. (48), which 
thus controls the low energy behavior of the FSC phase. 
Not only will phase-slips occur only in tightly bound 
pairs that will not be apparent at low energy scales, but 
phase slip dipoles will bind into quadrupoles, and isolated 
dipoles will not be evident at low energies. 

When both and the r]\ — the smallest dipoles that 
are least costly — are relevant, they will grow, inducing 
the other {?7s} to grow as well, until the small fugacity 
expansion breaks down. Below this scale, there will be 
free QPS's and no superconductivity: the resulting nor- 
mal phase is best studied from the weak coupling limit 
as we do in the next section, Sec. V. 

The behavior when r\\ is relevant but ^ is not is more 
subtle: either the SO* phase or the normal phase can ob- 
tain. If the phase-slips still bind into dipoles, but the 
dipoles do not bind into quadrupoles, the SC* phase 
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FIG. 7: Strong coupling phase diagram in K-Rq/R plane for various values of r. The three phase boundaries between the 
FSC, SC* , and NOR can be local, global or mixed. They are marked as following: (G) marks global (solid line), (L) - local 
(dashed-dotted line), and (M) marks mixed (dashed lines). The FSC-NOR transition is partially local (ML) and global (MG) 
as marked in the figure. The meeting point where the FSC-NOR boundary changes its nature is a bi-critical point marked by 
the gray dot. The open dot Marks the meeting of the three phase boundaries, and is a multi-critical point. Black dots mark 
where superconductivity breaks down in the limit of if — > 0. Note that the phase diagram in the vicinity of the sharp corners 
(bicritical points) may not be exactly correct due to large renormalizations of K which we otherwise neglect. These regions may 
give rise to quite complicated multi-critical behavior. The critical behavior in the vicinity of the bicritical and mult i- critical 
points is richer than near the other phase boundaries, and is characterized by very long crossovers. 
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obtains. In this phase we can ignore the ohniic dis- 
sipation and therefore also the dipoles as separate ob- 
jects: the proHferated dipoles essentially render the con- 
version resistance r infinite. Once dipoles proliferate, 
the SC* phase may still be unstable to single QPS. An 
instability to C once dipoles are proliferated drives the 
SC* phase to the NOR phase, but this stability must be 
considered differently. 

The qualitative discussion above gives rise to the main 
features of the phase diagram of Fig. 7: the exception 
being the weak-coupling section of the mixed transition 
that occurs for small K and small R\ this we discuss in 
the next section. Sec. V. 

The relevance of C, or ri\ about the FSC limit triggers 
the phase transitions from this phase to the normal and 
SC* phases as well as controlling the bicritical point at 
which the mixed phase transition changes character. In 
the limit of very strong coupling, the locations of the 
phase boundaries are given by the vanishing of the linear 
eigenvalues of C, and ryi. 

The SC* -normal phase boundary can be understood 
straightforwardly from the SC* phase (in which dipoles 
have proliferated) purely in terms of individual QPS. 
This transition is triggered by relevance of Q. Its loca- 
tion, however, cannot be precisely obtained in terms of 
the original parameters, even in the limit of strong cou- 
pling, because of the effects of the proliferation of dipoles. 



A. FSC- SC* phase transition [line (L)] 

In both the FSC and SC* phases single phase slips 
are irrelevant. Therefore, to determine the transition be- 
tween the FSC and the SC* phase wc need to know when 
■qi (the least costly dipole) changes from irrelevant to rel- 
evant about the Gaussian fixed line. The primary role of 
the non- linear RG flows for these purposes, is the gen- 
eration of {-qs} from C,"^ in Eq. (63) so that these will 
appear and must be considered even when they are ab- 
sent initially. [Note, however, that the scale at which 
rji becomes important is determined by the term so 
these second order terms can be crucial for understanding 
finite-temperature and other crossovers.] 

The first order terms of Eq. (63) for the flow of rji : 



(64) 



control the behavior near this transition. In the limit of 
infinite coupling they yield a boundary line between the 
FSC and the SC* phases: 



(65) 



which is thus a condition on the critical resistances, e.g. 

i?^(r) for this local transition. The quantitative super- 
fluid properties drop out as is the case for individual junc- 



tions. As wc shall see, for strong but finite coupling, the 
situation with a chain of junctions is more complicated. 

When r]i is irrelevant about the Gaussian fixed mani- 
fold, one must investigate whether there is a critical value 
of the fugacities above which rji grows. For /? close to 
unity, an expansion in the fugacities is still possible. The 
crucial terms are the creation of 772 by rji and the feed- 
back of r]2 on 771 . The structure is 

dm 
dl 



?7i (1 - /3) -h Cii2r?ir?2 



(66) 



^..(f-2«[f 



w 



^2117?! 



(67) 



with the Cs coefficients. Since the (linear) eigenvalue of 
ri2 is negative, there is a critical fixed point for [3 > \ at 



^J(3-l and 773 ~ /3 - 1 . 



(68) 



It can be seen that the higher moment dipoles have cor- 
responding fixed point values 

(/3 - ly'^ 



Vs 



(69) 



and can thus be neglected to leading order in (3—1. Since 
r]i will be of order the bare after the brief initial tran- 
sients, we see that for — R small but positive, there 
will be a critical value of C for the local transition 



a~(i??-i?)^ 



(70) 



using, for convenience, R as the control parameter; with 
r varying instead, the behavior is similar only with a 
different coefficient. 

Just above the critical fugacity, the crossover energy 
scale below which the chain will no longer appear critical 
goes to zero near the critical point as 



~ic-ar • 



(71) 



The critical exponent /U varies continuously with the re- 
sistive parameters. For /3 — 1 small, 

1 



Because this transition is essentially local — with dy- 
namic scaling exponent 2; = 00 — the diverging length 
scales can be subtle. But there should be some crossover 
in the spatial correlations at a length scale 



(73) 



beyond which the superconducting correlations behave 
somewhat differently. 

For this discussion, we have not been careful about 
where the RG flows go when the dipoles proliferate. If 
they go to the SC* phase, the above obtains. But it is 
also possible that they go near the SC* regime only to 
flow away at lower energies if C, is relevant. In this latter 
case, the transition will be from the FSC directly to the 
normal phase, as we discuss below. But we first consider 
the effects of dipole proliferation. 



21 



B. SC* phase and screening of dissipative 
interactions 

In the region of parameter space where rji is relevant, 
the simple perturbative RG flow, Eqs. (63), breaks down, 
and phase-slip dipoles proliferate. But this is not in itself 
sufficient to destroy superconductivity; the lead-to-lead 
coherence can persist as long as the quantum phase slips 
remain bound in pairs. [58, 83] Dipoles do, however, al- 
low the phase of a single grain to jump by 27r relative 
to its neighbors. This is the SC* in which each grain no 
longer has a well defined phase. But the global super- 
conductivity survives as long as C is irrelevant about the 
fixed line that controls the SC* phase. 

In the SC* phase, the dissipation no longer plays an 
essential role. This can be most easily seen by comparing 
the dipole term with the quadratic parts of the action for 
the V field in the Sine-Gordon action. At low energies, 
the fluctuations of the field 6 are small enough that one 
can ignore the 9j+i — Oj in the dipole term and approxi- 



Where we have neglected 6 — ijj cross-terms proportional 
to k^. These do not contribute to or alter the singu- 
lar interaction between phase slips in this regime. Note 
that the plasmon part of the interaction (the 9"^ part) 
has been renormalized by the phase-slip dipoles, but it 
is not screened, i.e., it still vanishes at u = k = 0. 
The dissipative part, however, does get screened reduc- 
ing the dissipative interaction between phase slips to 
Sdissix, t) oc 1/t^ when |t| > j;^. This should be 
compared with Eq. (42). Because of the finite screening 
time of the dissipative interaction between phase slips, 
this does not play a role for the asymptotic behavior of 
either the SC* phase, or the SC* to normal transition. 
The screening of the dissipative interaction by dipoles 
can be understood simply by circuit methods (Sec. II B) 
or by considering the phase-slips as a Coulomb gas. It 
plays a crucial role for finite JJ chains for which the ef- 
fects of dissipation are not completely eliminated. This 
will be discussed in a future publication [70] . 

If the system flows to the SC* fixed line when the dipole 
fugacity grows, then the critical behavior of the FSC- 



mate it by 

T]i cos[^/'j+i - tpj] . (74) 

When rji becomes of order O, the highest remaining fre- 
quency, the dipole term is comparable to the ii|a;|| V'j (tt') P 
and the r|a;||V'j+i(w) — tpj{uj)\'^ term at frequencies of or- 
der fi. [If r <C i? the comparison is somewhat more 
subtle.] 

Below this crossover frequency, the dipole term 
strongly suppresses fluctuations of neighboring ipj at low 
frequencies: together with the i?|cij||V-'j(i^)P term, it cuts 
off the divergent fluctuations of ^ while 9 continues to 
fluctuate similarly to in the absence of the dipoles. At 
this point, we can take the continuum limit without 
worry, replace the lattice with a high-momentum cutoff 
and expand out the dipole terms (analogous to treating 
screening in a metal as if by continuous charges — the 
Debye-Hiickel approximation) as 



(75) 



(76) 



SC* transition is that found in the previous subsection. 
Before returning to the trickier question of what happens 
when the flows eventually go away from the SC* regime, 
we briefly discuss the direct SC* -NOR transition. 

C. SC* — Normal transition [line (G)] 

Once dissipation has been eliminated, we are left with 
a simple low-energy effective action for the SC* phase: 
the standard chain of grains with self capacitance and 
connected by Josephson junctions. In the Sine-Gordon 
representation, we have 

•S'sc* ~ 

1 1 + I drC cos {9,{r)) (^7) 

where the tildes over the couplings are a reminder that 
they include renormalizations from integrating out the 
effects of the dipoles as in Eq. 76. 

The action Eq. (77) is most naturally studied by an 
isotropic RG rescaling both x and r, and integrating out 



r] cos a 



dij){x, t) ^ d9{x, t) 



« 

2 



dijj{x. 
dx 



d9(x, t) 
dx 



dx ' dx 
An approximate action for the SC* phase is then 

- E fdrC cos {ei{T)+ Mr)) 
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the modes of 9 that have large lo or large k. By carrying 
out the RG flow in x and r directions, the usual K-T RG 
equations for K = \/rJ and ^ obtain [84] : 



dK 
dl 

dC 



f = C(2-f). 

The SC* -normal phase boundary — 
— thus occurs for ^ <C 1 at 



(78) 



a global transition 



(79) 



When the initial is not small, the critical value of K will 

be increased both by the direct effects of ( and also via its 
generation of dipoles ({r/s}) and their later modification 
of K at intermediate scales as in Eq. (76). 

In the normal phase near this global transition, the 
characteristic time and length scales both diverge expo- 
nentially as exp[l/\/KG — K\ — and similarly in R—Rq 
when dissipation plays a role at intermediate scales and 
can control the transition. This contrasts with the power 
law divergences near the FSC-SC* local transition. 

Although we have here used the conventional isotropic 
RG to study tlw^ SC* - normal transition, it is important 
to understand whether the anisotropic RG we use in the 
rest of this paper can reproduce the KT results in the 
appropriate regime. We address this issue below in Sec. 
IV H. 



D. FSC — Normal transition [line (M)] 

We now turn to the most subtle transition: the mixed- 
character transition from the FSC to the normal phase. 
As the phase slip fugacity is increased, the superconduct- 
ing stiffness, K, is reduced, or the resistances increased, 
a quantum phase transition between the FSC phase and 
the normal phase of the chain can occur directly without 
an intermediate SC* phase. This mixed transition can 
be driven by two distinct mechanisms even in the strong 
coupling limit in which the Sine-Gordon description is 
useful. If single QPS's become relevant when dipoles are 
still irrelevant about the FSC fixed line, they will drive 
the transition inducing proliferation of both dipoles and 
isolated QPS as the fugacity C, grows. 

But if small dipoles become relevant when single QPS 
are still irrelevant about the FSC fixed manifold, it is 
nevertheless still possible for the dipolc-drivcn transition 
to be between the FSC and normal phases. This occurs 
if the flow towards the SC* phase from the proliferating 
dipoles is interrupted by the relevance of individual QPS 
about the SC* fixed line. We first discuss the latter case. 



E. Dipole driven mixed transition [line (ML)] 

The RG flows for this dipole driven transition are 
somewhat complicated. As shown earlier, if the QPS 
fugacity is irrelevant about the FSC fixed manifold when 
the dipole fugacities start growing — such as when the 
bare fugacity, C,, is above a critical value — C, will con- 
tinue to decrease until the effects of the dipoles are strong 
enough to change its flow. Crudely, this will happen when 
the rjsC, terms in Eq. (63) become large enough to dom- 
inate over the {\K + a — 1)Q term. As this will happen 
when some of the {ry^} are of order the energy scale fi, say 
at energy scale i?x , the system will soon after approach 
the SC* fixed line as discussed in Sec. IV A (other opera- 
tors neglected perturbatively will also become important 
on these scales). But by this point, if the system is close 
enough to critical, C, will have become extremely small. 
Thus even when it is relevant about the SC* fixed line 
and hence will turn around are grow on scales smaller 
than i^x , C ■^ill not become large enough to make indi- 
vidual QPS proliferate, until a much lower energy scale, 
-Eat < E-^ . 

Physical properties should exhibit both of the energy 
scales just beyond the FSC-normal transition: At high 
energies, the system will appear critical — FSC-SC* crit- 
ical — at intermediate energies, <^ E <^ E^, it will 
appear to be in the SC* phase — albeit with a phase stiff- 
ness that is too low to sustain superconductivity at long 
scales — and at asymptotically low energies, E <C En, 
normal behavior will obtain. This behavior is character- 
istic of a dangerously irrelevant operator — in this case 
C — about the critical fixed point. The energy scale E^ 
will go to zero as a power of the distance from criticality, 
but its exponent will be larger than that of Ex by an 
amount controlled by both the flow of ( towards the FSC 
fixed manifold, and its flow away from the SC* fixed line, 
as well as the flow of r/i away from the critical fixed point. 
The behavior of the length scales will be similarly com- 
plicated and the dynamic critical exponent relating En 
to the superconducting correlation length, will be larger 
than unity. 



F. QPS driven mixed transition [line (MG)] 

In the regime in which the FSC-normal transition is 
driven by QPS — the fundamental topological excita- 
tions — the phase boimdary and nature of the transition 
can be studied from Eqs. (63) for the flow of C,. For this 
transition, the dipoles play only a secondary role. 

In the region near the FSC fixed manifold, the pertur- 
bative RG equations (63) are valid and, in the regime of 
interest, the {rjs} are renormalized to zero and can be 
ignored. We thus need only the first two RG flow Eqs. 
(63). Working in units of the cutoff energy scale Q, and 
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length scale a, we have: 



dn 
dl 



dl 



-K, 



(80) 



a 



Equations (80) are similar to the Kosterlitz-Thouless (K- 
T) flow equations, except for the anisotropic scaling of k 

and J. Wc first proceed naively and transform to vari- 
ables in which the flows appear more isotropic, defining: 



J=Je-\ 
whereby Eqs. (80) become: 



(81) 



dl 



= 0, 



(82) 



dc 

dl 



C(i-iVKJ-a) 
Shifting and rescaling the variables via 



'-,/ 



(83) 



so that X parametrizes deviations from the special point, 
the flow Eqs. (82) assume the canonical K-T form 



dx 
dl 



^2 dz 
' dl 



-ZX. 



(84) 



The flows of Eq. (84) suggest a phase boundary be- 
tween the normal and FSC phases. In the FSC phase 
z flows to zero whereas in the normal phase z diverges 
and K flows to zero. The critical line separating the two 
phases is thus 



Xm 



(85) 



When C 1 the FSC-normal phase boundary in Eq. 
(85), X « 0, translates to 



K 



M 



3- 2a. 



(86) 



Although the phase boundary of Eq. (86) is correct for 

the renormalized K, for non-zero initial QPS fugacity 
C, the critical value of K will be slightly larger — by 
£)(C^) when the bare fugacity is small — because of the 
renormalizations of Eqs. (82). 

On the normal side of the transition, the behavior is 
subtle. The subtlety arises because the transformation 
between C and C, (Eq. 83 and 82 ) involves the scale of 
the RG exphcitly: C e~'/^C oc e~'/^2, so that Q may 



decrease even though z is formally relevant. But as z 
grows to the order unity, it will substantially decrease x, 
concomitantly, J, and hence the further growth rate of Q 
from Eq. (81). The growth rate of Q will then become 
rapid enough to overcome the e~'/^ factor and make C, 
grow as well. Thus after initially decreasing because of 
the factor e~'/^, C eventually turns aroimd and then be- 
comes large (at not much lower energies) at which point 
the behavior will be characteristic of the normal phase. 

The dipole fugacities, {r/g}, will be small as long as C, 
is, but when C, becomes large, they will grow (for large s 
as well as small). This implies that Josephson junctions 
are no longer independent of each other, and phase slips 
become spread-out over several junctions. In principle, at 
this stage we could switch to carrying out the RG isotrop- 
ically, but by then other higher order processes will also 
come in and expansions in powers of the fugacities break 
down. 

The critical behavior near this QPS driven mixed tran- 
sition will be similar to near a conventional Kosterlitz- 
Thouless transition, with, for example, energy scale, , 
vanishing exponentially in the inverse square root of the 
distance from criticality. The corresponding exponent 
is = oo in contrast to the finite but variable ^ that 
characterizes the dipole driven mixed transition. Associ- 
ated with the Ex will be one or more diverging length 
scales. There should be a scale that grows as l/E^ cor- 
responding to dynamical exponent z = 1, but there may 
be a second scale associated with the onset of substantial 
screening of k and spatially (rather than temporally) sep- 
arated QPS that will occur only after the dipoles come 
into play. We will not explore this issue further here. 

To make the above analysis solid. It is essential to 
show that the perturbative RG equations on which it 
is based are valid out to scales at which K{t) starts 
to decrease significantly from the screening by pairs of 
QPS. The rescaling from C to C was chosen to make these 
renormalizations, as in Eq. (82), small if and only if Q 
remains small. This could better (but equivalently to 
the needed order) have been done by defining instead 
C = ^fj C thereby making d log J/ d£ depend only on the 
combinations C and K = sfTk: the latter (in contrast to, 
e.g. c = \J J / k) is of order imity in the regime of inter- 
est near the transition. With this choice, it means that 
the dipole fugacities, {ry^}, are induced by {C,Y /J which 
decreases rapidly after the initial transients that create 
small dipoles. Thus the {r]s} should not be problematic 
as they will be decreasing and only give rise to relative 
changes: terms of order C^rjs . 

The primary potential problems are additional renor- 
malizations of J (or k) that are not small when Q is small. 
These can arise, for example, by the combination of a 
tripole — a QPS "dressed" by a nearby small dipole — 
and an opposite sign QPS. Such terms are equivalent to 
terms of order in other RG schemes, but these, and 
other terms higher order than quadratic, do not occur 
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directly with the RG scheme wc use: they are generated 
from other operators as in this case. This has the ad- 
vantage of making a cataloging of all operators and how 
they renormalize each other linearly and quadratically, 
equivalent to "all orders" arguments in other perturba- 
tive RGs. One can readily show that the generation of 
tripoles from a dipole and a QPS, and the subsequent 
renormalization of the tripole fugacities by dissipation 
plus the stiffness K, will result in effects of the tripoles 
back to C and to J that decay with energy scale while C, 
(and equivalently z) is decreasing. Further analysis along 
these lines shows that indeed, the result for the critical 
Km, Eq- (86) is correct. What is rather surprising, is 
that this arises from a basic (power-counting-like) renor- 
malization of C that is midway between the '1' for purely 
temporal rcscaling, and the '2' for space-time rescaling. 

By showing that the RG remains controlled out to 
this stage, we have verified that the FSC-normal phase 
boundary is in one regime indeed given by Eq. (86) in the 
limit of small fugacities; more generally the location of 
the transition will be exactly = 3 — 2a if the renormal- 
ized and hence measurable — low energy parameters, 
K and a, are used. 

Nevertheless, near the phase boundary on the normal 
side, the fact that C, first flows toward zero and only at a 
later stage becomes large, gives rise to a crossover in the 
resistance vs. temperature. From the last of Eqs. (80) 
we see that this will occur as long as 

K>2-2a. 

Thus the crossover occurs in the range of parameters 

2-2a<K <^-2a; (87) 

the physical consequences and origin of this region is dis- 
cussed in Sees. VII and II B. 



G. FSC — Normal bicritical point 

As shown above in Sees. IV E and IV F, the FSC- 
normal transition can have one of two distinct charac- 
ters: either a dipole driven transition with power-law 
singularities [line (ML) in Fig. 7], continuously variable 
exponents, and non-trivial relations between length and 
time scales, or driven by single QPS, with exponentially- 
rapidly decreasing energy scales and at least some space 
and time scales scaling similarly [line (MG)]. In the limit 
of small (bare) fugacities, the transition occurs where 
in the K,a,(i space the system is at intermediate en- 
ergy scales. The two critical surfaces in this space meet 
along a bicritical line, across which the critical behav- 
ior changes, although the phases on both sides are the 
same. At non-zero fugacity, C, the behavior is similar 
with the critical surfaces (actually now hypersurfaces) 



and location of the bicritical line (actually now a two- 
dimensional manifold) changed. When the critical fugac- 
ity is small, as occurs when /3 — 1 (see Eq. 64) is small 
(and positive), the behavior near the bicritical line can 
be treated perturbatively. If the multicritical manifold is 
approached along the dipole-driven critical manifold, the 
former can be found by considering the eigenvalue of C 
at the point on the critical fixed manifold, parametrized 
perturbatively by 77*. This eigenvalue will be increased 
by 0{rji)'^ from its value at the corresponding point on 
the FSC fixed manifold (all {-qs} zero). When the crit- 
ical eigenvalue, rather than 1 — — a passes through 

— ^, the critical fixed manifold will become unstable to 
individual QPS at low energies. This condition thus char- 
acterizes the multicritical manifold. For t]i just below the 
bicritical manifold the flow will be to the FSC phase, but 
because of the initial effects of rji on ^, will end up closer 
to critical than it would have been. Associated with the 
multicriticality will be complicated crossover behaviors 
that we will not explore. 

H. SC* — Normal transition revisited 

In the standard treatment of isotropic 2d, or 1+1, X-Y 
models, we use isotropic RG to obtain the K-T transition. 
This is also the way in which we analyzed the Global (G) 
NOR-SC"^ transition in Sec. IV C once dipolcs arc prolif- 
erated. As mentioned earlier, for consistency, we should 
be able to analyze this global SC* -normal transition by 
the anisotropic RG wc use in the rest of this paper. This 
also serves to make a more convincing case for the pecu- 
liar behavior of the QPS driven FSC-normal transition. 
Analyzing the global transition by the anisotropic RG is 
possible but clumsy; we only outline the basics here. 

It is useful to first consider what modification there 
would be to the small-fugacity flow equations if opposite- 
sign QPS on different junctions could annihilate each 
other under renormalization rather than combining into 
dipoles. Because of factors of = exp[— a/A] — sim- 
ilar to those in the generation of rjs by ("^ in Eq. (63) 

— only QPS pairs separated by less than or of order 
A{£) = J/k/H, will be renormalizcd away. But this is 
a factor of A/a more combinations than were included 
from QPS pairs on the same junction: the renormaliza- 
tion of J of order will thus be larger by a similar factor. 
We now recombine parameters similarly to what was in 
the mixed case: to make the renormalizations of J small 
when C is small and K of order unity. This now requires 
C = JC (contrast with ^/j ( in the mixed case) which will 
change the eigenvalue of ( to 2 — — a while keeping K 
marginal. In the absence of dissipation, this woifld yield 
a KT transition with critical value Kg = 4 as obtained 
from the isotropic RG of Sec. IV C. What is not imme- 
diately clear, is at what point the flows have crossed over 
into the normal regime and ^ cos 9j can be well approxi- 
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mated by C(l — ) and the Yl approximated by ^ / dx. 

3 

This can be checked by estimating the mean-square fluc- 
tuations of 9j using the resulting quadratic continuum 
form. As at low energy scales 1 /k is large and 1/ J small, 
this will be dominated by wave- vectors of order 1/A or 
smaller. The remaining frequency integral will be cut off 
at ^/QJ and thus dominated by frequencies of order f2 
when C = CJ/^ is of order unity as would be naively 
guessed from the rescaling that gave the KT flows. 

But this is not the whole story. In addition to the 
renormalizations of J, there should be similar (asymp- 
totically) renormalizations of k which controls the spa- 
tial gradients oi 0. In the anisotropic RG, these arise 
from replacing the dipole operators in the absence of ijj, 
cos{ej+, - 0j), by 1 - 2{ej+, - 6'j)2 « 1 - 2s^a^{d^e)^. 
This is justified if is sufficiently slowly varying in space. 
With the anisotropic cutoff, this is justified for s < A 
as, because of the low frequency cutoff relative to the 
wave- vector dependence, {{0j+s — 0j)^) ~ 2iir|s|/A in 
this regime. Thus the renormalizations of 1/k to order 

are of the correct form. By summing over the appro- 
priate range of s, their amplitude can be found to be of 
order Q'^h. which have the additional A factor that com- 
pensates exactly for the linear growth of 1/k with the 
inverse of the energy scale in the anisotropic RG. Sim- 
ilarly, the argument in the previous paragraph for how 
J is renormalized, can be justified by the replacement of 
{dr[0s+j + 0j\f cos(0j+s - 0j) by ^{dr0^s+jY renormal- 
izing 1/J. 

For the SC* - normal transition when the fugacities 

are small, the analysis above can only be done once the 
dipoles have been induced and proliferated enough to 
suppress the fluctuations of V and enable it to be ne- 
glected as in Sec. IV C . Again, the intermediate regime 
of the flows, where the {rjs} are of order unity, cannot be 
handled in a controlled manner. But this occurs over a 
small range of energy scales and thus will only result in 
factors of order unity. 



I. Multicritical point 

The three phase boundaries, FSC— SC* , SC* — 
normal, and FSC — normal come together at a multicrit- 
ical point as shown in the schematic phase diagram. Fig. 
7. The behavior near this multicritical point involves an 
interplay between the SC* -like behavior when dipoles 
proliferate, and the FSC-like behavior when the dipole 
fugacities are small. How these come together at the 
multicritical point is likely to involve flows in the inter- 
mediate coupling regime in which controlled calculations 
are beyond the scope of the methods of this paper. 



V. WEAK JOSEPHSON COUPLING LIMIT 

In this section we analyze the JJ chain when the 
Josephson interaction is weak: Ej <C Ec- In this limit 
the starting point is to assume there is no phase coherence 
between grains so that the currents between grains are 
predominantly normal. We can then consider Cooper- 
pair-tunneling events perturbatively. In this limit the 
Josephson coupling is sufficiently weak that superconduc- 
tivity can only be established if the dissipation is suffi- 
ciently strong to drive the system into the FSC phase: the 
SC* phase cannot occur as its stiffness K = 2'k^JEj/Ec 
is too small. We thus can explore only the dissipation- 
dominated transition directly from the normal to the 
FSC phase: as both local and global superconducting 
coherence arc established by this transition, it has mixed 
character. As discussed in Ref. 85, pair tunneling events 
are dual to quantum phase slips. Nevertheless, the treat- 
ment is substantially simpler than that of the mixed tran- 
sition in the strong coupling limit analyzed in the previ- 
ous section. 

A. Pair-tunnel events 

The basic objects that need to be considered are 
Cooper pair-tunneling events between grains. Their ac- 
tion can be formulated in the Coulomb-gas representa- 
tion. This allows the RG flows when pair-tunneling is 
rare to be worked out perturbatively. 

Consider first the Hamiltonian of a single capacitively 
shunted Josephson junction: 

Hjj = ^Q^-Ejcos(l>. (88) 

In the weak coupling limit, the charge fluctuations will 
be small relative to the phase fluctuations, and the junc- 
tion can be described in terms of charge states with wave 
functions: 

iS){(j)) = e'*"^, (89) 

where q is the charge imbalance across the jimction. The 
cosine term in the Hamiltonian thus induces hopping 
of Cooper-pairs, since e'"^ is a translation operator that 
changes the charge imbalance by 2e. We can describe 
the normal-to-superconductor phase transition of a resis- 
tively shunted Josephson junction in term of pair-tunnel 
events: when these are suppressed at low energies, the 
junction is well described by the charge state, Eq. (89), 
and is in the normal state. But when pair-tunnels pro- 
liferate, superconductivity is established and the phase 
becomes the good quantum number. In this dual de- 
scription, in opposition to that in terms of QPS, the 
quantum fluctuations are responsible for the supercon- 
ductivity rather than destroying it. 
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FIG. 8: (a) Cooper pair tunneling event, between two neighboring grains (solid line): the basic process for weak Josephson 
coupling. Charge relaxation is via the resistance 2r + R (dotted line), (b) Two nearly simultaneous pair-tunnels involving 
a common grain: the current through the conversion resistance r of that grain cancels giving rise to an effective interaction 
between the two pair-tunnels. These two pair tunnels can renormalize into a pair-tunnel over two junctions (c). 



To explore the pair-tunnel physics, we can start with 
the action in terms of the superconducting phases, {'/>j}, 
of the grains in the chain. This is given by Eq. (12). 
Since we are interested in the low energy dynamics on 



small length scales, we take the limit of small to of the 
quadratic parts of the action, keeping the spatial discrete- 
ness but dropping the 0{uj'^(f>'^) charging energy. The par- 
tition function in this limit can be written in the form: 



Z 



D[(f)iiT)]exp 



dk 



27r 27ra (i? -h r (2 - e}>"^ - e-'^"")) 



^iA(...)r+i:/..K,,».(A.M, 



(90) 



or equivalently with a local part on each grain, and an exponentially decaying — with w'^ — interaction of opposite 
sign between phases on grains s apart that exactly cancels for the total interaction (i.e., vanishes at fc = 0). This is 
because it is the phase differences across the junctions. 



Ai(T) = (/),;+i(t) - (t>i{T). 



(91) 



that enter here. Note that in the normal phase — the starting point here - it would be more proper to write the 
resistive part of the action in terms of d/S.j / dr which is proportional to the electrochemical potential differences across 
the junctions as the phase variables are not defined modulo 27r. The absence of superconducting aspects other than 
the cosine terms is because the quadratic action explicitly does not include the Josephson coupling. At low energies 
and when the capacitative interactions dominate, the densities of states parameters, and Z?s, drop out. Note that 
if there were no normal conduction across the junctions, i.e. i? = 0, the quadratic action would separate into a sum 
over non-interacting grains. 

The role of pair-tunneling events, and hence the superconducting behavior of the JJ chain, can be seen by carrying 
out an expansion of Eq. (90) in powers of Ej. The n'th order term corresponds to the probability weight of n 
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pair-tunnels occurring in spacc-timc. If the weight of the partition function (90) is concentrated in the term of O'th 
order in Ej, it means that no superconductivity is taking place at T = 0, since pair-tunnel events are suppressed. On 
the other hand, if the weight is concentrated on high order terms, pair- tunnel events proliferate — i.e., they happen 
frequently and the chain is superconducting. 

The expansion in Ej leads to a Coulomb gas representation of the action of pair-tunnel events and the partition 
function: 
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where {a = ±1} are the "charges" — i.e. signs — of the tunnels, the sum is over all distinguishable neutral 
configurations with ^ (t„ = 0, and the interaction between tunnel events is 



G{x, t) = -2qiq2 log 
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Rq 



Rq 



{5x,a + 5x,-a) 



(93) 



with 5x1, X2 the Kronecker 5. The short time cutoff, tq ~ 
{R + 2r)C/h is now the time duration of a phase slip, 
essentially the RC relaxation time of a charge imbalance 
between two neighboring grains. This and other cutoffs 
from the high energy physics, we approximate by a UV 
cutoff ilo — I/tq. The tunneling fugacity we denote 

6 cc Ej, (94) 

and it has units of frequency. 

The distance dependence of the interaction Eq. (93) 
is easy to understand. If two pair-tunnel events hap- 
pen through the same Josephson junction, the strength 
of the logarithmic interaction between them is 2-^^^ the 
total resistance that the Cooper-pair has to go through 
in order to relax back to an equilibrium charge distribu- 
tion (Fig. 8(a)) as found for a single shunted Josephson 
junction (see e.g. Ref. 58). If two pair- tunnels happen 
on neighboring jimctions, the Cooper-pair relaxation cur- 
rent overlaps on the common grain and thus involves the 
normal-to-superfluid resistance within the grain, r (Fig. 
8(b)). This reduces the logarithmic interaction strength 
of 2-^ between neighboring tunnel-events. Events sep- 
arated by more than one junction have independent re- 
laxations and hence no interaction as in Eq. (93). 



If two tunneling events that involve a common grain 
occur close together in time, the charge changes of the 
common grain can cancel and the combination is equiva- 
lent to a single pair-tunneling between further separated 
grains (see Fig. 8(c)). We denote by ^„ the fugacity of 
a pair-tunnel between grains separated by sa. Although 
in the original partition function, Eq. (92), ^„ = for 
n > 2, the long range pair tunnels will be produced dur- 
ing the RG flow. 

The strength of the mutual and self interaction of pair- 
tunnel events will determine whether they proliferate and 
induce superconductivity in the chain, or whether pair- 
tunnels form bound pairs and annihilate each other, keep- 
ing the chain normal. For wcak-Josephson-coupling, this 
can be analyzed via a perturbative RG analysis in powers 
of the {^J. 



B. RG Flow of the pair-tunnel fugacities 

To study the RG flows, it is easiest to go back to the 
initial action in terms of the superconducting phase dif- 
ferences, Aj, between the junctions: 
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As explained in the previous section, this is equivalent to 
a Coulomb gas of pair-tunnel events because 

exp I XI ) " ^« exp(0j+s - ^j) (96) 

\n=0 J 



increases the charge on grain j + s by 2e, and reduces 
the charge on grain j by 2e at rate ^g- We have explic- 
itly included longer range tunneling as this will, in any 



28 



case, be produced under renormalization from the basic 
nearest neighbor tunnehng which has rate = Ej. The 
dissipative interaction between a pair-tunnel of range s 
and a tunnel in the reverse direction a time r later, is 



2r + sR 
Rq 



log(r/To) 



(97) 



proportional to the total resistance through which normal 
relaxation must occur to compensate for a Cooper-pair 
tunneling between grains separated by s. 

The RG flows in powers of the fugacitics of pair-tunnels 
are obtained similarly to those in the strong coupling 
regime in terms of QPS (see Sec. VIIB of Chapter 5 of 
Ref. 80): indeed they have a similar structure to the RG 
for the various QPS dipoles in the absence of individual 
QPS (Sec. IV H). 

Including the combining of pair tunnel pairs - such 
as one from j — s to j together with one from j to j + m 
into a single one across m + s junctions — the RG flows 
are 
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Rq 
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(98) 



Note that we have neglected the renormalization of two 
overlapping pair-tunnel events which do not share any 
grains. Such events renormalize into a complicated non- 
contiguous compound four-tunnel events, and are more 
strongly suppressed. 



Phcise diagram and Normal-FSC transition for 
weak coupling 



flow is then toward large ^ and eventually to the FSC 
fixed manifold where the QPS fugacity, ( is zero. 

Because the non-linear contributions to the flows are 
positive, superconductivity can still occur when R + 2r > 
Rq and the {^s} are irrelevant about the Gaussian fixed 
point. If 



R + 2r 
Rq 



- 1 



(101) 



is small and positive, will be weakly irrelevant while 
the other {^s} arc still strongly irrelevant. Nevertheless 
(as in Sec. IV H) the feedback of these to important. 
For small this is dominated by the creation of ^2 at 
order £^f, and the feedback of this into d^i/d£ via the ^1^2 
term. There is thus a critical fixed point with 



(102) 



and ^2 ''^ S <^ ^i. It can be seen that the longer range 
{^s} attain fixed point values ^* ^ (5"/^. [Indeed, with the 
simple structure of the flows to quadratic order, an exact 
fixed point can be found: oc {By with B an r and R 
dependent factor that is of order \/5 for small S. But this 
is not a controlled expansion for the fixed point except 
when S is in any case small, or perhaps in some special 
part of parameter space.] Equation (102) corresponds 
to a critical fixed line that controls the critical surface 
of the mixed character normal — FSC transition. For 6 
small, the critical Ej ~ S.iiR,'''): when Ej < Ej 

the chain is in the NOR phase, while when Ej > Ej 
the chain is superconducting. Note that equivalently, one 
could pass from normal to superconducting phase at fixed 
Josephson coupling by decreasing R through its critical 
value Rm, which, in the weak coupling limit is 



The trivial fixed point of the RG equations for the pair- 
tunnel fugacities is the gaussian line with no pair-tunnel 
events: 



^s = Vs : 



(99) 



this controls the low energy behavior of the normal phase 
of the chain. The normal fixed point is stable as long as 
all the {^s} are irrelevant: since ^1 is the least irrelevant, 
this obtains when: > 1 , so that with R as the 

tuning parameter, the critical point in the limit of zero 
coupling is given by 



Rl,=RQ-2r. 



(100) 



When R < R^^, the normal fixed point becomes unstable: 

Cooper pair tunneling events proliferate and give rise to 
the superconducting phase. Since all the non-linear terms 
in the RG flows of the {r]s} are positive, these can only 
accelerate this flow so that the condition R < i?^^ is 
certainly sufficient for superconductivity Formally, the 
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(103) 



Thus the Josephson coupling decreases the amount of dis- 
sipation that is needed to cause superconductivity. This 
is the same underlying physics as that that makes the 
control parameter in the strong coupling limit ■^K + a 
which is the sum of a Josephson part, and a dissipa- 
tive part since a oc Xj \JR? + 2rR. But, not surprisingly 
given the basic electrical properties of the chain discussed 
in Sec. II B, the combination of parameters that enter 
in the strong and weak coupling limits are quite differ- 
ent. In terms of the superconducting stiffness parameter, 
K = 2iT^J Ej / Ec, the critical value of K for the mixed 
transition is 



Km ~ {R + 2r-RQ)^ 



(104) 



when this is small. 

The critical behavior near the FSC-NOR transition 
will be characterized by the exponent /i with the 
crossover energy scale going to zero a,s E^ ~ \Ej — Ej\'^ . 
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FIG. 9: Phase diagram with the weak Josephson coupling results. The solid lines represent the transition lines obtained with 
the strong coupling phase-slip approach (Sec. IV). The dotted line is obtained from extending Eq. (102) by solving the flow 

equations (98): K ~ --l=\/5- The two approaches seem to cross into each other continuously. The meeting point of the weak 
and strong-coupling boundaries is likely to be a bi-critical point where the transition changes its nature, but our methods are 
limited in that regime. 



The critical exponent, /i is a contimionsly varying func- 
tion of the resistances. For small 6, the above analysis 
yields 
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The rapid fall-off with distance, s, of the fixed point val- 
ues, (_* , of the pair tunneling amplitudes suggests that the 
fixed point does not have a long length scale associated 
with it, at least coming from the normal side. Whether 
there is a more subtle diverging length scale — e.g. that 
characterizes the decay with distance of the supercon- 
ducting correlations, we leave for future investigation. 

The important features of the weak-coupling phase di- 
agram are shown in Fig. 9. When 2r + R> Rq, the nor- 
mal phase can exist; otherwise the chain will be fully su- 
perconducting. But even for this high resistance regime, 
strong enough Josephson coupling can cause supercon- 
ductivity. Along the phase boundary — actually a crit- 
ical hypersurface in the parameters of our model — the 
critical behavior will vary continuously. 



This analysis of the weak coupling regime thus almost 
completes the phase diagrams of Fig. 9. The one excep- 
tion is the nature of the bicritical point at which the weak 
coupling section of the mixed-transition phase boundary 
meets the small-intermediate coupling KT section. As 
mentioned in the introduction, this does not appear to 
be amenable to perturbative analysis. A related question 
is whether, for intermediate coupling in terms of the QPS 
fugacities and pair tunneling amplitudes as well as K, the 
weak coupling segment of the mixed transition boundary 
can meet the large-intermediate coupling dipole driven 
segment without an intervening KT section. And, if so, 
do they come together smoothly or with another type 
of bicritical point. Again, this is not analyzable by our 
methods. 



VI. NANOWIRES 

The discussion so far has concentrated on the phase di- 
agram of a discrete chain of shunted Josephson junctions. 
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FIG. 10: Schematic of a QPS in a continuous wire. The 
size of the core of a QPS — in which the superconductivity 
is suppressed — is of order the superconducting coherence 
length ^. The wire can be considered as the limit of a chain 
of grains whose separation, a, is much smaller than f . 
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FIG. 11: Dipoles formed of a phase slip and anti phase slip 
that are a distance A ^ ^ apart. Their net effect is two 
partial phase slips of size A, a distance ^ apart. 



A. Dipoles in nanowires 



One of the motivations for this work, however, came from 
the observation of superconducting-normal transitions in 
nanowires reported in Rcfs. [23, 24, 86-88] (for ear her ex- 
periments see also [66-69]). In this section we will show 
how one can describe the superconductor-to-normal tran- 
sition in superconducting nanowires by carefuUy taking 
the continuum hmit a ^ of the discrete model dis- 
cussed before. Thermally activated phase shps in super- 
conducting wires have been considered by Langer and 
Ambegaokar [89] and by McCumber and Halperin [90]. 
Quantum phase slips in superconducting wires have been 
discussed previously in Refs. 64-66, 91-94. 

The fundamental difference between a continuous wire 
and a discrete JJ chain is the existence of the lattice con- 
stant in the chain. The interplay between the lattice con- 
stant, a, and the charge relaxation length, Aq, played an 
important role in JJ chains. In a wire, the characteristic 
dissipative length becomes 



Aq = ay/r/R = 



R/c 



(106) 



where 7 = 7^ is the phenomenological parameter de- 
scribing the conversion conductance-per-length from the 
normal to superconducting fluids, and p = R/a is the 
resistance-per-length of the wire. Nanowires are often 
modeled as discrete chains of Josephson junctions with a 
lattice constant of the order of the superconducting co- 
herence length, a = ^. [64] The motivation for this is 
that a typical size of a QPS should be of the order of ^. 
But from our discussion of JJ chains, it is clear that an- 
other important length is the smallest possible distance 
between QPS. In the case of JJ chains it is a, but for con- 
tinuous wires it may be arbitrarily small. So the correct 
description of the wire can be obtained from JJ models 
only by carrying out the a — ^ limit with the coherence 
length, fixed (see Fig 10). 



Physically, the continuity of the wire and the limit 
a ^ allows the formation of dipoles of phase slips that 
are arbitrarily close to each other — up to charge dis- 
creteness effects which are discussed below. The closer 
together the phase-slips are, the better is their screening 
from the dissipative interaction, which is the only low 
energy effect that suppresses the proliferation of dipoles 
in a discrete JJ chain. This decreased dissipative inter- 
action between them means that small dipoles in a wire 
will always proliferate. At zero temperature, a wire will 
thus always be either in the SC"^ phase or in the normal 
phase. 

This point can be made clearer by specifying how the 
continuum limit of a wire is taken. One method of taking 
this limit is spreading a phase slip over about ^/a ^ 1 
junctions, since a phase slip in a thin (diameter less than 
^) wire suppresses the superconducting order parameter 
over a length ^. This can be taken into account in the 
sine-Gordon action by changing the cosine term in Eq. 
(48): 



cos 



cos(V,(t) + 6',(r)) ^ 
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(107) 



where h{x) is a weighting function concentrated in a re- 
gion of length of order ^ with the constraint 



dxh{x) ~ a h[an) — 1 



(108) 



This constraint on h ensures that the total phase involved 
corresponds to one phase slip: it is forced by the global 
symmetry 0^9 + 27r. The UV cutoff will no longer be 
set by the shortest lengths, but be of order fi = c/^ or 
the Ginzburg-Landau time, which is related to the gap 
(in a way that depends on whether the superconductivity 
is in the dirty or clean limit). Although the operator in 
Eq. (107) breaks a phase slip into parts that are fractions 
of 27T, it does not lead to independent fractional phase- 
slips; because of the 27r periodicity of the cosine, the 
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fractional phase slips always appear together as a part 
of a complete phase slip. The structure of a smeared 
phase slip is shown in Fig. 10. The smearing means 
that the induced electrochemical potential drop caused 
by the phase slip is split into fractional potential drops 
that occur over roughly ^/a junctions adding up to the 
quantized drop of a single phase slip. The continuum 
limit a ^ can now be taken. 

Because of the dissipation, an individual QPS in the 
wire will have action that diverges logarithmically at low 
temperatures, as in the chain. The effective dissipation 
parameter is given by 



action 



a = J dx h{x) j dy h{y) 
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There are now two limits of interest. If Ag <C ^ (but still 
Aq ^ a so that the discreteness does not matter) the 
parameter is 



Rq 
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with the coefficient depending on the form of h{x) as 
should be expected: there is no unique way of defining 
the core size of a phase slip. This result involves the 
normal resistance over a distance ^: thus in the wire 
(rather than pa) plays the role of R in the JJ chain. But 
in general, is not the relevant resistance, in contrast 
to what would have been guessed from this analogy. 

If Aq ^ ^, ^ does not play a role in the dissipation 
and the behavior is similar to the JJ chain for Aq ^ a 
so that 



(111) 



analogous to l/2vri? for the chain. In this limit, a will 
be of order unity when the normal resistance of a length 
^ of the wire is still small. 

By analogy with the JJ chain, one would expect that 
the condition for breakdown of superconductivity would 
be controlled by the coefficient of the total interaction 
between QPS at the same position but different times, 
K + 2a. From the arguments of Sec. IV F, one would 
expect the critical value of this coefficient to be three. In- 
deed, there should be a change in the behavior near when 
this condition is met, but it will not be a true transition 
because of the role of dipoles. 

A pair of opposite sign phase slips separated by a dis- 
tance A — a dipole with moment A — will have infinite 
action at zero temperature because of the dissipation. 
But the dissipative interactions between the QPS that 
extend out to A of order the greater of Aq and ^ will 
reduce the coefficient, /3(A) of the logarithmic action of 
the dipole. 

As in the chain, a pair of equal and opposite dipoles 
with the same location but separated in time will have 



2/3(A)ln(T/ro) 
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with the dissipative parameter a function of the resis- 
tances and the dipole moment, A, relative to ^ and Aq. 
In the limit that A ^ min(^, Aq), for smooth weighting 
function h, the action can be expanded in powers of A 
yielding 
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Rq 
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For ^ <C Aq, there will a change in behavior from 
quadratic to linear dependence on A for A ^ ^: the 
larger moment behavior then being like the JJ chain with 
Aq <C a so that (3 saturates to « 2q; at A ~ Aq. For 
^ P will crossover from quadratic to saturating at 
A~ ^. 

Because (3 decreases to zero as A 0, independent of 
the dissipative parameters or of ^, for sufficiently small 
A, P will be less than unity and thus dipoles will prolif- 
erate — as long as very small dipoles can be modeled as 
we have done here. The characteristic moment A below 
which dipole proliferation occurs is 

Ax ~max(AQ,Ov/pe/i?Q (114) 

provided this is less than ^, otherwise a different depen- 
dence obtains. Nevertheless, we see that independent of 
the relative magnitudes of ^ and Aq, dipoles of moment 
^ will not proliferate as long as l/y^pTr < ChRq with 
Ch an order-unity coefficient that depends on details of 
the structure on scale ^. 

If dipoles of size of order ^ can proliferate, the behavior 
will be similar to the J J chain when s = 1 dipoles. the 
smallest size, proliferate. The dipole proliferation will, 
at low energies, lead to screening of the dissipative inter- 
actions between individual QPS and thus to either the 
SC* phase or the normal phase depending on whether 
the renormalized K is greater than or less than four. 

But if only dipoles of size A <C ^ can proliferate on 
their own, the behavior will be quantitatively different. 
Although it is not clear what pairs of QPS whose sep- 
aration is much less than ^ represent, they should cor- 
respond to large amplitude localized fluctuations of the 
superconductivity that act to decouple the superconduct- 
ing from the normal degrees of freedom. The effects of 
these will be to renormalize the dissipative interactions 
between larger moment dipoles and individual QPS by 
decreasing the rate of the processes, parametrized by 7, 
that relax imbalances between the superconducting and 
normal voltages. Again, if the dipole moments can be 
arbitrarily small, this will eventually lead to either the 
SC* or the normal phase. But if the bare rate of the 
small-dipole-like processes, t?(A), is small, there should 
be interesting crossovers as temperature is lowered as- 
sociated with the growth at low energies of these dipole 
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fugacities before they fully screen the dissipative interac- 
tions. 

In principle, even in a chain of grains linked by Joseph- 
son junctions QPSs need not take place only across the 
junctions but could happen inside a grain — albeit with 
much higher action and hence exponentially lower fugac- 
ity. But in the limit of extremely low temperatures, such 
processes need to be considered and the chain will be 
more like a nanowire. In particular, the effects of QPS 
dipoles with very small size within one grain will, by 
the arguments above, be relevant. Eventually, these will 
screen the superconducting fluctuations from the dissi- 
pative interactions and drive the system into the SC* or 
normal phases: again, if arbitrarily small dipoles can oc- 
cur, the FSC phase will not exist in the limit of asymp- 
totically low temperatures. 

B. Minimum size QPS dipoles? 

In all the above, we have assumed that arbitrarily small 
QPS dipoles can exist: in the continuum model we have 
used, even with short distance cutoffs, there seems to 
be no reason these cannot occur. But as the QPS are 
quantum objects, we must be careful of the possibility of 
complex — rather than purely real — actions. 

In a QPS dipole, the superconducting phase, (l){x), 
inside the dipole winds by 27r as r goes from — oo to 
+00 (see Sec. IIIB). As long as (j){x,T) is periodic in 
re (0, a QPS dipole of the opposite sign is needed to 
satisfy this condition. But in reality, the true condition 
is that (j) is periodic in to /3 modulo 2n. In the limit 
of zero temperature, this will have no effect in the action 
we have studied. But in general, there is an additional 
term in the action, analogous to a Berry phase, of the 
form: 

Syj = i J dx J drhs^ , (115) 

with fis having units of number density. In a Galilean 
invariant system, at zero temperature fis would be the 
density of Cooper pairs, equal to half the electron number 
density. But more generally, it can have any value, pos- 
itive, negative, or zero: ns is some effective density for 
the superfluid that has no known simple interpretation. 
It will depend on particle-hole asymmetry and other fac- 
tors that are not understood [95, 96]. For vortices in a 2d 
film, fig corresponds to a dual magnetic field, and it gives 
rise to a vortex hall effect [97]. [Note that in Ginzburg- 
Landau theory, only the combination ng/m* is physical: 
ns, in contrast to the phase stiffness, T often called 
ps — has no physical meaning. Quantum mechanically, 
this is no longer true.] 

The action of a dipole of moment A has an imaginary 
part of its action from the Berry phase term: 

g,Mpoic ^ 27rmsA . (116) 



This suggests that dipoles for which n^A is not an inte- 
ger will undergo destructive interference implying that A 
should be quantized in units of 1/ns- The QPS dipoles 
can be considered as a vortex which first crosses the wire 
in one dircetion, and then returns and cuts the wire in 
the opposite direction at a distance A away from the first 
crossing, before returning to its origin. Interpreting hs 
as Cooper pair density, the quantization of A means that 
the vortex had to "go around" integer numbers of Cooper 
pairs [98]. 

If there is indeed a minimum size QPS dipole with 

A„i„ = J- , (117) 
ns 

and the low frequency coupling to the "normal" degrees 
of freedom is still effectively dissipative for such small 
dipoles, then it is possible that the dissipation could be 
sufficiently strong to suppress these smallest dipoles. [Al- 
though, the effective coupling, (3{Amin), would be deter- 
mined by the shorter length scale physics which will not 
be well parametrized by 7 and p.] If the smallest dipoles 
are indeed suppressed, then an FSC phase — and a tran- 
sition from it to an SC* or a normal phase — would be 
possible in nanowires. 

It is worth noting that effects of a small local distur- 
bance can be long range in time and give rise to actions 
that diverge logarithmically as temperature is lowered. 
An example is an impurity in a metal that is moved by 
a small amount: the divergent action is associated with 
the "orthogonality catastrophe" caused by the changes 
in the electron wave-functions induced by the altered po- 
tential. [99] Better known, but closely related, is the 
X-ray edge singularity [100]. Thus it is not unreasonable 
to think that a local change in the superconducting de- 
grees of freedom could result in a logarithmically infinite 
action associated with the response of the quasiparticles 
and other low energy degrees of freedom to the change. 

Whether or not the FSC phase can occur, even in prin- 
ciple, we leave as an intriguing open question. 



VII. RESISTIVITY 

In the previous sections we concentrated on the zero 

temperature phase diagram of the RSJ.I chain shown in 
Fig. 1(c). We now extend this analysis to discuss the 
temperature dependence of the resistivity of the system 
using the RG flow equations obtained in Section HID. 
One of the most surprising results is that in the normal 
phase close to the superconductor-to-normal transition 
the system exhibits a non-monotonic R(T) dependence. 
In this regime one flnds that the resistivity flrst decreases 
with decreasing temperature, then saturates and stays 
nearly constant at a value R{T) <C R, i.e., much smaller 
than the normal resistance of the chain, down to very 
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low temperatures, and only then finally starts to in- 
crease. During this long crossover, the Josephson junc- 
tions themselves form a resistive channel, parallel to the 
shunt resistors. Since experiments have a lower limit of 
temperature, such behavior may appear to support the 
existence of the "metallic phase" [14] separating the su- 
perconducting and insulating ones. 

We start by showing how one can obtain the resistance 
vs. temperature curves from the RG flow for the QPS fu- 
gacity. We then proceed to discuss the finite temperature 
resistivity in the FSC and SC* regions, and explain the 
origin of the quasi-metallic behavior in some of the NOR 
region. 




A. Resistance of a single junction 

The RG analysis presented in this paper allows one 
to calculate the scaling behavior — including crossovers 
— of the temperature dependence of the resistance of a 
JJ chain. The simplest way to do this is to carry out 
the RG flow until the UV cutoff becomes of order the 
temperature. At lower energy scales the renormalized 
phase slips are only weakly interacting and their effects 
independent. 

To obtain physical quantities, we must distinguish be- 
tween the "bare" parameters — those that appear in 
the action with the initial high energy cutoff and the 
renormalized parameters. In situations in which there 
might be confusion, we denote the renormalized param- 
eters with subscripts, i, corresponding to energy scale 

We first consider a single resistively shunted Joseph- 
son junction and estimate the (superconducting) voltage 
fluctuations across it at a low temperature, T. At low 
frequencies, the resistance of the superconducting part 
controls the low frequency fluctuations via 

{V{oj)V{uj')) = RT ■ 2n6{oj + oj') . (118) 

The voltage fluctuations caused by QPS are tricky to 

estimate because of the interactions between them. But 
if we renormalize until an energy scale at which the QPS 
are weakly interacting, T{£) ~ Cl, then the fluctuations 
can be estimated simply. 

The phase-slip fugacity Q and the temperature of a 
single RSJJ obey the flow equations [85]: 



b. 
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dl 
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with the high energy cutoff, fio, initially of order the 
plasma frequency of the circuit. When the RG flow 
reaches the stage in which T; ~ O, only one pair of 
phase slips (they must come in pairs because of the pe- 
riodic boundary conditions in r) is likely in the rescaled 
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FIG. 12: Semi-log plot of the resistance vs. temperature 
for various values of K. (a) Resistivity in the SC* phase. 
K ^ 3, 3.2, 3.4, 3.6, 3.8, 4, 4.2. (b) Resistivity in the FSC 
phase. K + 2a = 2, 2.25, 2.5, 2.75, 3, 3.25, 3.5. In both plots 
the resistivity at the bare critical value of the stiffness appears 
as a thick black line {K = 4, and K + 2a — 3). The non-zero 
phase-slip fugacity shifts these lines slightly from criticality. 
Resistivities in the non-mouotonic pseudo-metallic region arc 
shown as thick gray lines. Temperature is normalized by the 
plasma frequency of an individual junction. The contribu- 
tion of the parallel channel of normal electrons to the total 
conductivity has been subtracted. We use the approximation 
~ g-^K/rr initial condition, which corresponds to the 
bare action of the QPS So = i-^EjjlEc [36, 101]. 



imaginary time duration l/T^ and these are, in any case, 
weakly interacting. At this scale, the renormalized fu- 
gacity is 



Ct = = ln[^^o/T]) . 



(120) 



The phase slip pairs will then occur with probability of 
order Ct/^^ time of order UjTt ~ fi. A phase slip 
contributes to ^(a; = 0) — a quantity that is scale inde- 
pendent — an amount h/2e. The rescaling of T and of 
5{ui + uj') in Eq. (118) cancel, and the renormalization of 
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C to Ct then yields 



2-2a 



(121) 



The resistance increases with decreasing temperature in 
the normal phase which obtains when Rq/R < 1, and 
goes to zero at zero temperature in the; superconduct- 
ing phase which obtains when Rq/R > 1. This re- 
sult, Eq. (121), is valid only for T < Hq. For higher 
temperatures, T > fig activated scaling of the form 
Rj.T ^ RQCxj>i-CjEj/T) with Cj = 0(1) obtains. 
How those match together can be understood in terms 
of the fugacity (.. This is of order fJo cxp(so/?i) with sq 
the action of a phase slip sq ~ Ejtq ~ Ej/Qq since the 
time scale of a QPS is tq. If the imaginary time h/T is 
smaller than tq, the action of the phase slip will be re- 
duced by a factor of order H/Tto- This gives rise to the 
activated resistance in terms of Ej/T. 

The measured resistance of the shunted junction is 



Rm^ ~ (Rjj) ^ + R 



(122) 



If in Eq. (121), C 



no 



or larger, the Josephson 



27tRq 

junction becomes essentially insulating, and most of the 
current flows through the shunt resistor, R. 

Note that at any non-zero current, some Cooper-pair 
tunneling will occur and non-linear resistance V/I, will 
no longer be given by Eq. (122). But these corrections 
vanish in the limit of zero-temperature and zero-current. 



B. Resistivity in the FSC region 

We now turn to the JJ chain. In the FSC phase, each 
junction behaves similarly to a single junction in its su- 
perconducting phase and thus contributes roughly inde- 
pendently to the total resistance. Phase slip dipoles, with 
fugacities {r/s}, are the explicit manifestation of inter- 
action between phase slips on different junctions. Thus 
when all the {rjs} are irrelevant about the FSC fixed line, 
the phase slips in each junction will be almost indepen- 
dent although their dissipation involves overlapping 
resistors. Thus each Josephson junction will act like the 
a single junction with resistance given by Eq. (122), with 
the appropriate the rcnormalizcd fugacity at scale 

T ~ r2 — of individual QPS in the chain; 



R 



JJ 



{2ey 



(123) 



At intermediate temperatures (what we mean by inter- 
mediate is discussed in Sec. IV F), Eq. (123) predicts a 
measured resistance per unit length, 



Rjj /a ■ 



(124) 



with the (positive) exponent F depending on the normal 
resistances and the superconducting stiffness, K. 



In the FSC phase, the exponent F is given by minus 
twice the eigenvalue C [from Eqs. (80)] 



T = K + 2a-2 



(125) 



Since the boundary between the FSC and normal phase 
when the {r/s} are irrelevant — the mixed transition — 
is given by Km = 3 — 2a, we see that on the critical line, 

in contrast to the single junction critical point at which 

the measure resistance will be roiighly temperature inde- 
pendent. The origin of the factor of T can be seen from 
the transformation that related the flow of Q to the flow 
of the parameters z Q cxp(— i?/2) and a; = + a — | 
which obey the Kosterlitz Thouless fiows Eqs. (84). 
When Q is small, a; = — at which z is weakly scale 
dependent — marks the FSC-NOR transition. The scale 
dependence of the relationship between z and Q trans- 
forms Eq. (123) to Rjj ^ RqTqC,^ ^ Rq^z^ giving rise 
to the factor of T in terms the "natural" variable z. As 
we shall see, this really is natural for the KT-like global 
transition from the SC* phase to the normal phase. 

The behavior of Rjj{T) near the FSC-NOR transition 
is shown in Fig. 12(b). The results in Eqs. (80-82) lead 
to a regime of parameters in which a crossover between 
quasi-superconducting and normal behavior takes place: 
in particular, the resistance exhibits a minimum at a low 
temperature. This behavior arises on the normal side of 
the transition, but still close to it: when 



. K 3 
^<T+"<2' 



(127) 



so that Q initially decreases with scale even though z is 
increasing. If the bare fugacity is small, this gives rise 
to an intermediate temperature regime with a power law 
decreasing resistance with exponent < F < 1. But at 
low enough energies, the rapid increase of z will give rise 
to proliferation of dipoles and individual QPS and the Q 
will start to grow. The junction then crosses over to in- 
sulating behavior as the temperature is lowered further. 
Near the transition to the FSC phase, the crossover tem- 
perature, Tmini at which the resistive minimum occurs 
becomes very low: 

Train ~ f^o exp {-h/^/KM - (128) 

with h a coefficient proportional to the bare QPS fugacity, 
when this is small. Below the crossover temperature, 
the measured resistivity, pm, of the chain will increase as 
the temperature is lowered, until it reaches the resistance 
of the shunt resistors, R/a. 



C. Resistivity in the SC* region 

The low temperature behavior in the SC* phase can 
be understood similarly. When the renormalized dipole 
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fugacities arc large, dissipation will be screened and we 
need only consider individual QPS and the interactions 
between them mediated by the superconducting degrees 
of freedom. 

It is now simplest to use the isotropic RG and rescale 
space and time. When the RG flow reaches the stage 
where T{£) = fl, the spatial size of a QPS is flo/T times 
its original size. Therefore, when resistances per unit 
length are found from renormalized quantities, a factor 
of the length rcscaling is needed to obtain the physical 
resistivities. The resistivity of the Josephson junctions, 
will thus be 



pjj - ^^icrm' - ^mr ( ^ ) (129) 



fin a 



T 



with 



r = 1 + (Jf - 4) 



(130) 



the K — 4: being minus twice the eigenvalue of ( in the 
isotropic RG. We have assumed that the basic length 
scale that enters is the grain spacing, o, and used the 
bare QPS fugacity (: more generally these will depend 
on the higher energy processes as discussed below. 

The resistivity as a function of temperature is plotted 
for several values of K in Fig. 12(a) Since for K < Kg = 
4 the chain will be normal at low temperatures and the 
resistivity increase to R/a, for the intermediate range 



3<K kA 



(131) 



the resistivity will decrease with dec;rcasing T and have 
a minimum at a crossover temperature Tmin given by 
the energy scale at which the QPS proliferate enough to 
screen the longer-range interactions between them driv- 
ing the junctions normal. This range is analogous to the 
intermediate regime of Eq. (127) near the FSC phase: In 
both cases the resistivity seems to show signs of super- 
conductivity before exhibiting insulating behavior. For 
ii' w 3 the flat resistivity curves over a substantial tem- 
perature range could be confused for a quasi-metallic 
phase. 



D. Roles of dipoles 

In the JJ chain, if the QPS fugacities arc small and 
the dissipative interactions substantial, the QPS fugac- 
ity can be irrelevant, but the dipoles may be relevant 
about the FSC fixed line. In this case, the SC* phase 
behavior can obtain only at low enough energies below 
which the QPS dipoles proliferate: this makes the tem- 
perature dependence of the resistivity much richer than 
that discussed so far. If the SC* phase is not stable, the 
resistivity near the dipole driven FSC-NOR transition, 
will involve a crossover from dipole to individual QPS 
dominated resistivity. 



The temperature behavior of the resistivity in both 
these cases splits into two regimes. Since the {rjs} are 
zero initially, and get generated at order at tempera- 
tures in the range To < T < CIq wc expect the JJ chain 
to behave as though it were in the FSC region. The 
crossover temperature Tp, away from the FSC regime 
is the energy scale at which 771 becomes of order il: at 
lower energies, the dissipative interactions get screened 
(see discussion above Eq. (75)). At energies higher than 
To, the scaling is similar to that of individual junctions, 
while at lower temperatures, it will become isotropic as 
in the SC* regime discussed above, and exhibit normal 
behavior at sufficiently low energies if K is too small to 
stabilize the SC* phase. 

We now estimate Tp, and discuss its consequences for 
the low temperature behavior. The flows of 771 when it 
is relevant can be well approximated by an initial rapid 
regime in which r]i becomes of order (^/Qq, followed by 
a scaling regime in which it grows with exponent 1 — /3 
with P = 2a(l — w) as in Sec. IV H. Thus it will become 
of order Cl at an energy scale 



no 

c 



2/(1-/3) 



(132) 



for (3 below its critical value of unity for either the FSC- 

SC* or the FSC-NOR dipole-drivcn transitions. This is 
valid as long as Q grows less rapidly than rji: i.e. if 
2-2a - K < 1-/3. This is always the case if K > 1 
and sometimes also at smaller ii' if w is substantial. If 
C| grows more rapidly than exp[(l — P)£], then r]i follows 
C| and individual QPS will drive the crossover to normal 
behavior. 

For T > Td, the individual-junction controlled behav- 
ior found above for the resistivity in the FSC phase (and 
for crossover away from FSC) will obtain: 



2a-2 



(133) 



But below Td the temperature dependence of the resis- 
tivity will be like that in the SC* phase (or the flow away 
from that) with p„ ~ T^'^. 

As before, if ii' < 4, at low enough temperatures the 
chain will crossover to normal behavior and the resistivity 
saturate at R/a, while if if > 4 the resistivity will go to 
zero more rapidly than linearly. 

Wc have seen in this section that the interplay between 
local and long-length scale physics gives rise to interest- 
ing behavior of the resistivities. Perhaps surprisingly, 
the temperature dependence near the critical point of the 
QPS driven mixed transition, for which the critical value 
of the parameter ^k-I-q; is |, midway between the values 
of 1 and 2 expected for local and global transitions, nev- 
ertheless is similar to that of the global transition with 
Pm ^ T, with the T arising from spatial rescaling. 

Unfortunately, the form of the various crossovers of 
the resistivity as a function of temperature, which would 
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FIG. 13: Possible realization of a two-fluid JJ chain: meso- 
scopic superconducting grains deposited on top of a metallic 
nanowire and separated by a thin insulating layer, to produce 
an SIS junction. 

yield firmer testable predictions, are not readily acces- 
sible to our analysis since these involve analyzing a 
sine-Gordon model in the intermediate coupling-strength 
regime. In addition, we have been rather cavalier about 
which resistivity we are considering: in general, even in 
simple scaling regimes, the resistivity at w ~ T and the 
DC resistivity can be quite different. The latter certainly 
requires proper real-time analysis. 

VIII. DISCUSSION 

In this final section we discuss some implications of 
the results presented in this paper, and issues associated 
with these. 

Chains of Mesoscopic Josephson Junctions. Sev- 
eral groups have recently reported experimental realiza- 
tions of mesoscopic superconducting grains [102] and ar- 
rays of low-capacitance Josephson junctions [19, 35, 36, 
47-51]. The possibility of introducing ohmic dissipation 
using shunt resistors has been demonstrated for a single 
quantum Josephson junction by Pentilla et al. [36], and 
Watanabe and Haviland [19], and for arrays of Josephson 
junctions by Takahide ei aZ. [53] and Miyazaki et aZ. [57]. 
A potential realization of a chain of mesoscopic Joseph- 
son junctions with ohmic dissipation is shown in Fig. 13. 
It consists of equally spaced small superconducting grains 
on top of a metallic nanowire with thin insulating layers 
separating the two materials. This system is a realiza- 
tion of the model that we presented in Sec. Ill, with 
the possibility of controlling separately the parameters 
Ej, Ec, R, and r by changing the width of the metal- 
lic wire, the size and separation of the superconducting 
grains, and the strength of the tunnel barriers between 
the grains and between the grains and normal wire. A 
schematic phase diagram of this system is presented in 
Figs. 7, and 9 (we did not consider odd-even effects in the 
grains, which must be irrelevant if r < Rq, as discussed 
below) . Both the phase diagram and the character of the 



transitions are affected by the superconducting-normal 
relaxation parametrized by r. [124]. 

If the metallic wire is highly resistive and/or the tun- 
nel barrier between the S (superconducting) and N (nor- 
mal) parts weak, the dissipation will not be important, 
the FSC phase will not occur, and the superconductor- 
to-normal transition will be determined primarily by the 
competition between Josephson coupling and Coulomb 
charging energy of the grains. This will be a quan- 
tum Kosterlitz-Thouless transition, isotropic in space 
and time, with the energy scale decreasing exponentially 
as the transition is approached. 

If the wire is less resistive and/or the SN tunnel junc- 
tion stronger, the FSC phase can occur. The FSC to nor- 
mal phase boundary and the FSC to SC"^ phase boundary 
will be dominated by dissipative effects. These change, 
in particular, the nature of the superconductor-to-normal 
transition which will have a mixed character with both 
local and global physics involved. This mixed transition 
can either be KT-like — although modified from the con- 
ventional KT transition — or have strongly anisotropic 
power law scaling with continuously variable exponents. 

The renormalization group analysis presented in this 
paper allows us to obtain the scaling form of the chain's 
resistivity as a function of temperature in the vicinity 
of the various superconductor-to-normal transitions; the 
results are presented in Fig. 12. Near the SC"^ to normal 
phase boundary, our results agree with the results in Refs. 
[64, 65], but the analysis of R{T) in the vicinity of the 
FSC-to-normal phase boundary is new. We find that in 
both cases the system has extended crossover regions on 
the normal side of the transition, in which superconduct- 
ing tendencies give rise to a weakly temperature depen- 
dent resistivity that decreases as a small power of T over 
a wide temperature range before turning up at the lowest 
temperatures as it becomes asymptotically normal. Near 
the superconductor-to-normal phase boundary this up- 
turn in R{T) occurs at temperatures that are are either 
power-law or exponentially small in the inverse of the de- 
viation from the zero-temperature transition, depending 
on the regime. 

Our RG approach to the problem is also new, and al- 
lowed us to expose a myriad of crossover effects occur- 
ring near phase boundaries. Particularly interesting is 
the interplay between global and local mechanisms for 
the breakdown of superconductivity. The mixed FSC - 
normal transition, for example, can arise from the local 
resistive environment, or from a combination of the long- 
wavelength plasmon dissipation and the local resistive 
parts. The phase boundaries for which these two mecha- 
nisms obtain meet at a bicritical point, which has a rich 
crossover behavior yet to be analyzed in full. The SC* - 
normal transition is shown to be global, i.e., determined 
by the long-wavelength excitations, and thus, for an in- 
finite chain, the low energy critical behavior is almost 
independent of the resistive shunting or the two-fluid na- 
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FIG. 14: Schematic of a type of superconducting nanowire 
with dissipative degrees of freedom: a superconducting layer 
is deposited on top of a metallic nanowire with an insulating 
layer in between. 



ture of the grains, except for their effects in modifying 
the effective capacitances and Josephson couplings. 

Other aspects of the phases and phase transitions could 
also be investigated for a chain of grains. Of particular 
interest would be tunneling into pairs of grains, studying 
the behavior as a function of their separation and of tem- 
perature, and how these change near phase transitions. 

Superconducting Nanowires. Various types of su- 
perconducting nanowires can have dissipation down to 
zero temperature. The simplest is perhaps a wire with 
regions that remain normal, so that it has a positive den- 
sity of states at zero energy. Another possibility is a 
sandwich of superconducting and metallic nanowires sep- 
arated by a thin insulating barrier (see Fig. 14). But even 
nanowires made of a conventional gapped superconduc- 
tor may have low energy degrees of freedom that give rise 
to dissipation when they are close to a superconductor- 
to-normal transition. The conventional BCS theory of 
superconductivity, which leads to an energy gap, is intrin- 
sically based on a mean-field analysis and thus unlikely 
to be reliable in the vicinity of a quantum phase transi- 
tion when virtual quantum phase slips occur. A roughly 
analogous effect is seen in superfluid ''He in a disordered 
porous medium (e.g. Vycor): an apparent normal fluid 
component has been observed down to the lowest tem- 
peratures [103]. In general, one would expect dissipa- 
tion of some sort — but not necessarily ohmic — near 
superconductor-to- normal quantum transitions as long as 
some of the Cooper-pairs break apart near the transi- 
tion. This may not occur if the Cooper pairs form tightly 
bound bosons that still exist in the non-superconducting 
phase, but more generally should occur. Furthermore, 
if the normal phase is a Bose or Fermi glass, there will 
be a constant density of states of low energy excitations 
and concomitant dissipation arising from local two level 
systems. That said, the approach we have taken in this 
paper is phenomenological: we assume that dissipation is 
present in at least some superconducting nanowires and 
then study its consequences. 

For superconducting nanowires, our analysis of the 
simplest model suggests that nanowires will always be ei- 
ther in the SC"^ or the normal phase with the FSC phase 



unstable to small QPS dipoles. The superconductor- 
normal transition would then generally be of global KT 
character. But if the fugacity of phase slips is low, 
and the dissipative effects large, there will be interesting 
crossovers in the temperature and length-scale dependent 
properties in both the superconducting phase, and near 
the transition. 

Nevertheless, as discussed in the Sec. VI, these con- 
clusions may not be correct. In particular, charge dis- 
creteness effects can lead to interference between QPS 
dipoles which would set a minimum size for dipoles in 
the limit of zero temperature. If the dissipation is strong 
enough — a condition that will depend on short-length 
scale physics in addition to the mesoscopic resistive pa- 
rameters — then it is possible that the FSC could be 
stabilized. These effects clearly need further exploration. 

Finite length wires. In the SC"^ phase at zero tem- 
perature, the only low energy effects of dissipation in our 
nanowire model arise from the global fc = mode and 
thus depend solely on the total resistance of the wire. In 
finite wires this would mean that at very low tempera- 
tures, when the wire acts like a single Josephson junction, 
its total resistance will determine whether it is supercon- 
ducting or resistive. At higher temperatures, T > c/L, 
the behavior of the wire could exhibit a K-T crossover 
which is tuned by the stiffness, or thickness, of the wire. 

We note, however, that in finite wires the resistivity 
may determine other parameters that are important for 
the superconductor-to-normal transition, such as the fu- 
gacity of the QPS [24] and the strength of quantum fluc- 
tuations [93]. So the superconductor-to- normal transi- 
tion in flnite wires is determined by both the resistance 
and the resistivity. A more detailed discussion will be 
given elsewhere [70]. We should note that other effects 
may also be important for the case of superconducting 
nanowires, such as the suppressed Tc due to Coulomb 
interactions as in Refs. 104, 105. 

The picture arising from our analysis conforms with 
and enhances that in Ref. 94. There it is assumed that a 
finite nanowire is shunted externally by a resistor. This 
resistor gives rise to a zero-temperature SC-NOR tran- 
sition, but the plasmon degrees of freedom in the wire 
give rise to sharp crossovers. Our picture puts the pic- 
ture of Biichler et al. [94] on firm microscopic footing by 
showing that, quite generically, the resistivity of the wire 
needs to be taken into account as a single resistor that 
shunts the entire length of the wire. 

Experiments on nanowires. The hypothesis that 
dissipation plays an important role in superconduct- 
ing nanowires is supported by the recent experiments 
of Bezryadin et al. [23, 24, 86], in which the 
superconductor-to-normal transition was observed to de- 
pend on the normal state resistance of the wire. Later 
experiments [24] on longer wires showed a transition that 
depended on the resistivity of the wire, rather than its 
total resistance. Bollinger et al. [87] recently carried 
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an extensive study of shorter wires 200nm), and 
observed results consistent with both previous experi- 
ments. At this point it is still not clear whether the zero- 
tcmpcraturc SC-mctal transition in finite-length amor- 
phous nanowires depends on their total resistance, or 
primarily on the resistivity. 

Charge discreteness effects. Before turning to 
broader implications, we briefly raise various issues. The 
most serious caveat to our results arise from aspects of 
charge discreteness that wc have ignored. In the formal- 
ism we use, the charge discreteness not only gives the 
basic periodicity of the superconducting phase, 4>, of the 
Joscphson coupling, but also implies periodic boundary 
conditions modulo 27r in the imaginary time path integral 
on (0, e.g., in Eq. (12). We have ignored the effects of 
non-trivial windings of the phase and thus the linear ^ 
term that only matters when (j){T + (i) ^ (I>{t). These can 
cause Berry-phase-like interference effects among QPS 
that wc have ignored. 

For single Josephson junctions, an approach to dissi- 
pation that respects charge quantization is to explicitly 
analyze quasiparticlc tunneling as carried out, for exam- 
ple, in Ref. 37. It is found that the quasiparticlc tun- 
neling has very similar effect to Ohmic dissipation when 
the (continuum) electrostatic equilibrium corresponds to 
an integer number of Cooper pairs on each electrode. 
The generalization of this to JJ chains in the absence 
of dissipation has been analyzed in Ref. 106 and others 
[107]: the simple model we use is only directly applica- 
ble when there is an integer number of Cooper pairs per 
grain (or other commensurate filling for which the grains 
can act in groups). At incommensurate filling factors 
the Cooper pairs are always delocalized in the absence of 
randomness. Hence, naively, one would expect that if the 
ohmic shunt resistors are replaced by quasiparticlc tun- 
neling and motion of discrete electrons, and the conver- 
sion resistance by quasiparticlc creation and annihilation 
processes caused by the dynamics of the superconducting 
phase, (j), the conclusions of our paper would only be valid 
for integer number of Cooper pairs per grain. This addi- 
tional effect of charge discreteness, which would severely 
restrict the applicability of our results, is most likely to 
be problematic when the shunt and conversion resistances 
are high: this will increase phase fluctuations and result, 
in any case, in charging eflfects destroying the supercon- 
ductivity as wc have found. When the resistances are 
low — the regime in which the FSC phase is predicted 

— it is less clear what the effects of charge discreteness 
are as the number of Cooper pairs on a grain is no longer 
a good quantum number. We leave for future work this 
issue, as well as the question of whether near to, but on 
the normal side of transitions, there is an intermediate 

— or fuller — temperature range over which our results 
apply. 

The case of effectively commensurate filling, corre- 
sponding to the coefficient of the d/dr, ns = is the 



case in which a transition is possible even in the absence 
of dissipation. [106] From the superconducting side, this 
should be driven by proliferation of individual QPS as 
the absence of resistive interactions, dipolcs have finite 
action, and thus will always occur as fiuctuations. 

Randomness. For another issue, the effects of ran- 
dorrmess, the above discussion is also relevant. Giamarchi 
and Schulz analyzed a strictly one-dimensional boson sys- 
tem with a random potential. [59] They discussed the 
important role of the random Berry's phase, which con- 
trols the equilibrium local density of bosons, and found a 
transition between the superfluid and the insulating Bose 
glass phases that is driven by pinning of the density fluc- 
tuations of the superconducting phase by the random 
potential. How wires with spatially random properties 
would behave in the presence of "normal" carriers and 
dissipation is unclear. A similar question arises for the 
KT-like transition predicted by Altman et al. [107] for 
a commensurately filled chain of grain with random self 
capacitance and nearest neighbor Josephson interactions. 

In the strong randomness limit of a JJ chain with dis- 
sipation, the superconducting transition is effectively lo- 
cal as we have found here, but with subtle interactions 
between the local phase slips. With strong randomness, 
the transition should thus be dominated by weak links on 
many length and energy scales, and may be controlled by 
an infinite randomness quantum critical point as found 
in other random quantum systems. [108-111] 

The effects of randomness on normal carriers need also 
be considered. Localization effects in a nanowire typi- 
cally become important below the Thoulcss temperature 
which is roughly the spacing of energy levels in a seg- 
ment of length such that its total resistance is of order 
Rq. Simple estimates of this for a J J chain suggest this 
temperature is of the same order as T* up to factors of 
R/Rq — which make the Thouless temperature much 
lower when this factor is small, i.e., strong dissipation. 
But with a system like that of Fig. 13, the tunneling bar- 
riers between the grains and the normal nanowire both 
play roles: in this situation, the conversion resistance, r 
will be large, and the Thouless temperature is likely to be 
well below T*. Near transitions to superconductivity, the 
superconducting fluctuations will affect the normal elec- 
tron transport and any possible localization. As these 
decrease the overall resistance, it is likely that they also 
lower the Thouless temperature — perhaps all the way 
to zero at the transition. This argument is supported 
by the results of Giamarchi and Schulz: the SC* phase 
is stable to randomness if its stiffness parameter, K, is 
sufficiently large. These, and other issues, surely merit 
exploring: the interplay between superconductivity and 
randomness is still largely an open field. 

Higher Dimensions. Finally, we briefiy discuss the 
implications of our results for thin films and bulk ma- 
terials. The percolation picture of the superconduct- 
ing to normal transition in two and three dimensional 



39 



systems of grains, [10, 39] suggests that close to the 
transition the conductivity wiU be determined primar- 
ily by quasi-one-dimensional percolation paths. This, 
and the local effects of dissipation on each junction, sug- 
gest that many of the effects that we study here for one- 
dimensional systems should apply to transitions in gran- 
ular systems more generally. These effects may lead to 
various interesting crossovers, including the possibility of 
non-monotonic temperature dependence of the resistiv- 
ity (see also Ref. [58]) and quasi- metallic behavior near 
to putative superconducting-insulator transitions as ob- 
served in a variety of systems. [14, 112, 113] One of the 
most intriguing (and poorly understood) features of such 
transitions in thin films is the existence of a "supermetal- 
lic phase" characterized by a small but apparently weakly 
temperature resistivity down to very low temperatures — 
extrapolating, it appears, to zero-temperature. [62] This 
behavior is in striking contrast with the theoretical pic- 
ture of the superconductor to insulator transition in two 
dimensional systems presented in Refs. [81, 106, 114- 
117]. In this, the resistivity in the limit of low temper- 



atures should go to either zero or infinity. The origin of 
the observed metallic behavior in thin films is still un- 
clear, although several groups have recently addressed 
this problem [16, 61, 118-120]. A possible origin of such 
a 'supermetallic phase' in thin films are the percolation 
aspects of the transition combined with effects analogous 
to those that we have discussed in this paper. In particu- 
lar, these should be significant when, near the transition, 
there are many low energy excitations that give rise to 
dissipation. 

We note that after this paper was complete, Ref. 121, 
which analyzes a similar problem with similar methods 
appeared online. 
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APPENDIX A: DERIVATION OF LOW ENERGY ACTION AND JUSTIFICATION OF 

APPROXIMATIONS 

In this appendix we outline derivations of various low-energy forms of the action that are used in the text, from the 
basic action Eq. (12) of the chain in terms of the superconducting phase and the electrochemical potential of each 
grain: 

Schain = Sq + S^^^ + S*^^ + Sj 

i 

+(3j:(j^\VN{u>n)-VscM')k 

w„ ^' ' ^ (Al) 

-!^dTY.{Ejcos{4>i+^{T)-Ur))) ; 

i 

at zero temperature, the sums over Matsubara frequencies can be replaced by integrals: j3 ^ / which we will 
use herein. 

From Eq. (Al), one can integrate out the normal voltages, Vn{x, r), to obtain: 

SoHa^n = 5, + i/ |^/f([c- F{k,^)\ \Vsc{k, u:)\' + ^^|^ \Vsc{k, i^f) (A2) 



where 



(A3) 
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with 



_ C{Dn+Ds) ^ 
*^ - Dn+Ds+CDnDs ~ ^ 

B - ^s*^ 



Dn+Ds+CDnDs ■ 

In expression (A2), the last term and the C term (which is a shghtly modified capacitance in the hmit CDm,s ^ 1) 
comprise the usual action for a dissipatively shunted JJ. Together with these two terms, the third, comphcated looking, 
term, with coefficient F{k,u)), gives rise to the two fluid behavior. It is from the frequency dependence of F that the 

basic energy scale 

arises. The low energy regime of interest is w <C T*. Note that because of the local nature of some of the important 
physics, we must take the low frequency limit at fixed k; if the limit is taken in the opposite order, the results can be 
very misleading. 

At low energies, the action becomes 

which, on approximating C by C, is the form of the action we use in both the strong and weak coupling limits. 

In the weak coupling limit, we drop the capacitative term as this is unimportant at low frequencies. 

In the strong coupling limit, the capacitative term is important and the Villain approximation for the Joscphson 
coupling and ensuing transformations can be carried out straightforwardly as discussed in Sec. Ill A. In general, this 
is only appropriate for low energies, in particular for u) <^ Qlr, the inductive relaxation rate 

2e^' 



flLR = y-^j EjR, (A7) 

since Ej is inversely proportional to the "kinetic inductance", Lj = (|^)^ ;g- of the junction. Note that ^lr is the 
same order as the plasma frequency of the junctions, reduced from this by a iactor of KR/Rq which is of order unity 
in most of the regimes of interest. Thus the Villain approximation is valid wherever we have used it. 

The action in terms of phase slips, the interactions between them, and the equivalent sine-Gordon representation in 
terms of the dual fields 6 and tf} can be derived straightforwardly. Note that we could have gone more directly to the 
sine-Gordon representation by decoupling the inter-grain terms in Eq. (Al) with the fields "0 and 9, then integrating 
out Vn and (f). The periodicity of (f) then gives rise to an integer constraint on {9 + tp) /27r which after integrating 
out the high energy fluctuations, becomes the cos(^ -|- V') in the sine-Gordon action. The appropriate low frequency 
approximations appear naturally in the intermediate representation. In particular, the quadratic coupling between 6 
and %l) that appears if this route is followed, is unimportant at low frequencies. 
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